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ABSTRACT 


Three-dimensional finite element analyses were 
performed to investigate the stresses and settlements during 
construction of Mica Dam. This high earthfill structure 
in a skewed steep valley was a suitable subject for studies 
to develop further confidence in simplified analytical 
approaches which can successfully model complex field 
situations. 

Despite the economically necessary coarseness of 
the analyses, representative behaviour matching field 
measurements and tests was calculated using a simple linear 
elastic approach. Conclusions are drawn as to the situations 
where such a simplified representation of soil behaviour is 
Valid. 

Confirmation of previous investigations was made, 
suggesting that the stress state in the core of the dam is 
very much less than expected. Some implications of this 
important feature are discussed, 

A series of parametric studies of a homogeneous 
Mica Dam showed that a Young's modulus of 7800 Ksf. anda 
Poisson's ratio of 0.28 yielded good estimates of shell 
settlements. Similar values for core settlements were 4500 
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Analyses of the zoned embankment were made using 
Simple linear elastic parameters, including bedrock settle- 
ments, and also a 2D plane strain condition was imposed. 

A simple approach to nonlinear behaviour was made by allowing 
deformation moduli to vary with stress level. Average 
valuesof Young's modulus and Poisson's ratio were 7800 Ksf 
and 0.28 for the shell, and about 2400 Ksf and 0.33 for the 
core, to match field settlements. 

It was demonstrated that the dominant deformational 
behaviour of Mica Dam during construction could be 
satisfactorily represented by 2D plane strain analyses. 
Stress distributions for various analyses are given to 


emphasize the low stress state in the core. 
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CHAPTER I 
INTRODUCTION: SCOPE OF THIS STUDY 
1.1 Role of Deformation Studies in Design and 
Monitoring of Embankments 
As part of the design process the engineer must 
estimate the movements, strains, and pore water 
pressures which may develop in his dam at various 
times in its life. He makes these estimates largely 
on the basis of laboratory tests and theoretical 
analyses, but since the theories are still ina 
state of development he can only rely on them to 
the extent that they have been demonstrated to check 
actual measurements on completed structures. 
Sherard, et al. Earth and Earth-Rock Dams, 1963. 

Once estimates are made it is the engineers' 
responsibility to ensure that the structure exhibits the 
behaviour predicted. If deviations from predictions are 
significant the questions "why" must be answered, particularly 
in the context of safety of the structure. 

Significant advances in analytical power have 
become available with finite element techniques, particularly 
in dealing with stresses, deformations, seepage, and 
consolidation. The wealth of information potentially available 
at the design stage of a project is staggering when the whole 
rather than the parts can be observed, namely with 3D 
analyses. The engineer's perception of his ability to know 


"what is going on" is being expanded greatly and personal 


adjustments need to be made to this new era. 
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Instrumentation programs have been widely used for 
some time now, but it is only in the age of the finite 
element where the full potential of behaviour monitoring 
can really be appreciated. The new analytical power is 
yielding far more than before so there is even more reason 
to check field behaviour to verify analyses. The require- 
ments and responsibilities demanded by society's new and 
growing environmental consciousness accentuate the need to 
know in detail "what is going on" with a project and make 
this knowledge freely available. This thesis aims at 
being a contribution to such need. 

In this chapter background concepts for the thesis 
are discussed. Succeeding chapters amplify the details 
and present the results, and conclusions are made specifically 
regarding this study on Mica Dam but more generally about 


the role this work is designed to fulfil. 


at2 Importance of Documented Case Histories 


The economy of earth and rockfill structures lies 
in the utilization of readily available materials, and 
their performance is a function of the engineering efforts 
PUT Rene gelectee design, and construction. Usually one 
Cannot transfersdatal from: sitesto» site. since each project is 
essentially unique. An analysis manipulates data from one 
site and in the process many assumptions, usually whole 


series of assumptions, are made. 
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3 
To verify that the analysis is valid measurements are 

made at the site, and the agreements made and questions 
raised remain attached to that site. To transfer these 
experiences of results, measurements, comparisons and 
questions it is necessary to document the whole situation 
thoroughly. Only by doing this can confidence.be born in 
the versatility of the analytical procedures. This process 
is accentuated by a continuing stream of theoretical and 
analytical developments. At the level where most engineers 
work numerical results are required. When an analysis produces 
these results, and documentation is there to support it, the 


general well-being of the discipline is advanced. 


1.3 Problems of Arching and Stress Transfer 

The overburden stress at a point in a soil mass is 
equal to the weight of materials in a vertical column of 
unit area above the point. It is independent of any 
stress-strain parameters in an ideal, semi-infinite 


continuum, and is also the maximum principal stress, 


O4 ee 2 ¥ hh 
where Dest maximum principal stress at a point 
Y = material unit weight 
h = vertical extent of the material above the 


point in consideration 
To obtain the other two principal stresses involves 


stress-strain theory, and if linear elasticity is assumed 
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where Ky = = (a special case really) 
U = Poisson's Ratio 
Jor53 = intermediate, minor principal stresses 
respectively. 


If the stress distribution in a soil mass is altered 
from this simple state by an imposed boundary condition, 
then in general the major principal stress will not equal 
the overburden stress and in particular it will usually be 
less than overburden. The difference may be considered as 
taken up in stress redistribution to satisfy equilibrium at 
the boundaries. The phenomenon of stress changes to satisfy 
equilibrium and compatability under boundary influences is 
loosely termed stress transfer, load transfer, or arching 
(since the principles are embodied in structural arches). 

In embankment engineering it is common to talk of 
"cross valley arching" since the process so closely parallels 
the distribution of stresses in an arch bridge across a 
valley. Furthermore if the embankment consists of two 
different materials and a softer material "hangs up" on 
a stronger one major principal stresses in the two decrease 
and increase respectively from the stress state present if 


they had identical properties, and we speak of "load 
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transfer"rors "stress transfer” .p* Since! it ws* useful! to 
separate the effects if possible, this thesis will 
consistently use the term "arching" to mean cross-valley 
stress transfer due to geometrical configuration of the 
boundaries of an embankment across a valley. "Stress 
transfer" will refer to the effects of variations in 
deformation behaviour, and will also be used where the 
geometrical and deformational effects cannot be separated. 
In three dimensions the stress state at a point is 
quantitatively represented by the octahedral normal and 


shear stresses, 
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If no stress transfer effects are present and all stresses 
are calculated using simple elastic theory (with only two 


independent parameters), then the same elastic theory yields 
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and any deviations from this ratio thus indicate the presence 
of boundary effects or differences in deformational behaviour. 
This is an interesting idea but becomes rather impractical, 
consequently any evaluation of stress transfer or arching 


is made in this thesis by considering the relative values 
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of overburden stress and major principal stress. In most 
situations where stresses are due to material self-weight 
only, the major principal stress is close in magnitude 


and often in direction’ to the vertical stress. 


1.4 Review of Some Stress Transfer Effects in Soils 

In practically all engineering situations stress 
transfer effects will be present. It is necessary to have 
a feel for the processes and effects since stress transfer 
can have important implications. 

Terzaghi (1943) devotes discussion to the simplest 
example of stress transfer, arching in this context. He 
clearly develops the concept that shear stresses are 
mobilized in order to satisfy equilibrium in a mass when 
boundary conditions are modified. Similarly, concepts 
useful in the understanding and design of braced excavations 
can be developed. (Terzaghi and Peck, 1967; Lambe and 
Whitman, 1969). Refer to Fig. 1.1. 

Several examples of arching and stress transfer are 
cited by Tschebotarioff (1973) in the context of consolidating 
fills, vibrations, and flexible bulkheads. He indicates 
clearly how the mobilization of shear stress in response 
to boundary conditions influences the overall stress 
distribution on the structures discussed, 

A theoretical study of the development of stresses 


over flexible and rigid culverts within a soil mass, including 
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Modifying boundaries (after Terzaghi, 1943) 
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Braced excavations (after Lambe and Whitman, 1969) 
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the effects of different construction methods, was reported 
by Costes .(1956). Centrak-toehis- analysis isthe .satisfying 
of equilibrium between different portions of the soil mass/ 
culvert system by development of shear stresses along 
boundary surfaces. Some rather arbitrary estimates of 
arching effects are also involved, and he justifies these 
on the basis of field observations. 

Stress transfer is always present in some form in 
embankments, particularly zoned fills. Kulhawy and Duncan 
(1970) performed plane strain finite element analyses of 
Oroville Dam, and indicate that the analysis explains quite 
well the construction behaviour of the structure. The stress 
transfer between core and stiffer shells, and around the 
concrete core block, were clearly shown. By plotting the 
shear stresses in terms of strength mobilized toward failure 
(using the Mohr-Coulomb criterion) they demonstrated the 
relatively high proportions of strength mobilized at the 
sites of stress transfer processes. 

Lefebvre, Duncan and Wilson (1973) compared 2D and 3D 
analyses of an idealized homogeneous earth dam, to illustrate 
some features of finite element analysis of embankments. 
Arching is clearly shown by the effects of valley slope 
steepness (Fig. 1.2). In a similar study Palmerton and 
Lefebvre (1972) demonstrated the strength mobilization 
against the valley walls accompanying arching, and similarly 


demonstrated stress transfer between shells and core of 
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Major Principal Stress Contours, Longitudinal Section 
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Maximum Shear Stress Contours, Longitudinal Section 


FIGURE 1.2 EFFECT OF VALLEY STEEPNESS ON ARCHING 


(after Lefebvre, Duncan, and Wilson, 197 3) 


tov baatbetsaro! eave dD mans laeloaizs vote 


SS 
if ee | ee 


nO — — 


La 
frase Teolbdakgr)..,wueme) nascad tenl® nannies 


as anche ae 


SOLSAP (PLANE STRAIN ANALYSIS) 


10 


SOLSAP (3-D ANALYSIS) 


50 


Tima 4 


r Tfailure x100 % 


FIGURE 1.3 STRENGTH MOBILIZATION DUE TO CORE-SHELL STRESS 
TRANSFER (after Palmerton and Lefebvre, 1972) 
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FIGURE 1.4 HYDRAULIC FRACTURE MECHANISM 


(after Kjaernsli and Torblaa, 1968) 
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diGierentypropertiess(FRigsel1.3)2 

Because of the stress transfer and arching effects, 
the stresses in the core of an earth dam are usually less 
than overburden. Schober (1970) presented a rigorous study 
of the stress distributions in Gepatsch Dam, where field- 
measured earth pressures in the core were at least 50% less 
than overburden, and relatively high shear stresses were 
developed at the core-shell interfaces. 

Squier (1967, 1970) discussed in considerable detail 
the instrumentation results from El Infiernillo Dam and 
Netzahualcoyotl Dam in Mexico. Various modes of load transfer 
were mentioned. During construction of both these dams the 
cores tended to settle with respect to the upstream and 
downstream shells, giving rise to certain patterns of 
movements and stresses. When the reservoirs filled and wetted 
the upstream shells, these settled with respect to the cores 
causing a reversal of the shear stress mobilization at the 
upstream interfaces and correspondingly complex movements 
in the structures. Stress concentration due to non-uniform 
abutment slope was also apparent. 

In what has come to be the classic example of the 
consequences of stress transfer in dams, Kjaernsli and 
Torblaa (1968) reported the development of horizontal cracks 
through the core of Hyttejuvet Dam causing leakage and 
requiring a remedial grouting program. Stresses in the core 


were so low that the hydrostatic pressure of reservoir water 
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13 
caused hydraulic fracturing. Although it may be quite 
difficult to visualize the mechanism by which the fractures 
occurredsthermfact is that’ they) did* occur) /and”theWauthors 
presentedParpossible explanationy(F2g?eo1e4). 

Stress transfer and arching is a ubiquitous feature 
of soil behaviour. The finite element technique is potentially 
capable of shedding light on the mechanisms involved and the 


way stress transfer develops over the history of a project. 


1-5 Status of Finite Element Analysis in Geotechnical 
Engineering 


Much time and energy has already been expanded by 
numerous authors in assembling state-of-the-art reviews of 
this subject. The finest examples available to this writer's 
knowledge are the series of addresses to the W.E.S. 
Symposium on Applications of the Finite Element Method in 
Geotechnical Engineering held at Vicksburg, Mississippi, in 
May 1972. It is easily obvious from the literature now 
available on the subject that quite enough techniques 
are available, but it is still necessary to prove the 
usefulness of the finite element as a design tool. 
Practical situations can often be solved using the simplest 
methods which will yeild reliable results of appropriate 
accuracy, if the analysis is based upon a realistic 
description of soil behaviour. 

Few authors have reported successful usage of the 


method during design of a project. In only a very few of 
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these were the analyses carried out prior to or early in 
the construction period. So one good question which 
practising engineers rightfully ask is "just how were the 
analyses performed?" It is always much easier to calculate 
the required answers after an event has occurred. The need 
for finite element research aimed at developing useful design 
tools is to predict behaviour before it happens. It is to 
be expected that the method will not be so successful as it 
appears now. How will consolidation effects be represented? 
How can representative deformation parameters be selected 
taking into account the effects of compaction machinery? 
Will simple (elastic?) parameters, carefully chosen using 
reasonable assumptions, give good enough answers or will 
sophisticated models using great arrays of test data be 
ecu ced? These are only some of the questions, and it is 
emphasized that they are only poorly answered after the event. 
Prior to application in design work some experience has to 
be accumulated in the analysis of case histories. The 
principle of performing "post mortems" is a natural 
development in any field where new concepts are being 
implemented in practice. 

Far from being content that we now have a powerful 
new design tool, we should devote as much effort as possible 
to documented case histories where the analysis has been 
done prior to construction. Skermer (1974) indicated that 


1969 analyses performed for Mica Dam produced construction 
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ahs, 
movements larger than were measured at the site. The most 
important areas of development of the finite element 
method in geotechnical engineering must be in improving 


our predictive capabilities. 


1.6 The Finite Element As a Predictive Tool in Design 

The finite element method can be used to compute 
pore pressure due to construction, consolidation, and 
seepage under complex boundary conditions. It can also be 
used to determine stresses and displacements at various stages 
of construction, under complex boundary conditions, and with 
simple or esoteric stress-deformation behaviour of many 
materials. 

At the design stage where configurations and economics 
are being evaluated many analyses may be required. Computer 
budgets may be relatively restrained so only very simplified 
analyses could be undertaken. There is considerable argument 
for one branch of research effort to develop influence factors 
for the combination of materials and configurations usually 
considered in design, as suggested by Covarrubias (1969). When 
a design is being finalized special analyses will be required, 
but many simplifying assumptions would still need to be made 
particularly regarding stress-strain behaviour of the materials. 

We are not yet at the stage where most simplifi- 
cations can easily be justified. There is a great need for 


further research to explain just what effect certain 
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simplifications will have. There is, above all; a great 
need for the persons who use the method in design to present 
their analyses and results, good and bad, in order to further 
confidence in the capabilities and limitations of the method 
in future design work. Otherwise there is the danger of 
finite element analyses being taken too seriously without 
good reason, Or not seriously enough, again without good 
reason. 

Significant contributions to the development of 
the finite element method to investigate cracking of earth 
dams have been made by Eisenstein, Krishnayya, and Morgenstern 
(1972, a,b) and Krishnayya (1973a). The 3D finite element 
program developed by Krishnayya (1973b) has been adapted by 
the writer to the analysis of construction of Mica Dam. The 
size and height of this structure, its location in a skewed 
steep river valley, and the materials used in construction all 
Make it an excellent case to study and it is hoped the results 
of these analyses will contribute to the understanding of such 


a situation in future engineering works. 
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CHAPTERS I 


INTRODUCTION TO THE MICA DAM PROJECT 


Pee iGeneral Description of MicasDam 


Mica Dam is a high earthfill structure sited ina 
narrow gorge where the Columbia River leaves the Rocky 
Mountain Trench about 80 miles north of Revelstoke, British 
Columbia. The reservoir created will extend for over 130 
miles in the Trench, from Golden to Valemount. 12 million 
acre-feet of live storage will be provided under the terms 
of the Columbia River Treaty (1964), and eventually 2500 
MW of electrical power will be generated for Canadian 
consumption. A general description of the project has been 
given by Meidal and Webster (1973). 

The damsite is in rugged, glaciated, mountainous 
terrain. Bedrock is a Precambrian metamorphic series of 
mica shist, granitic gneiss, pegmatite and white marble 
(in decreasing order of occurrence). Isoclinal folds are 
present, with bedding and foliation planes generally dipping 
downstream (southwards) at about 20°. Pleistocene geology 
is represented by glacial deposits and volcanic ash layers, 
river terraces and slope deposits. The sound bedrock 
conditions would be suitable for various types of dam, but 
after costs and material performances were evaluated a zoned gravel 


and rockfill structure with a near-vertical central impervious 
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18 
glacial till core was chosen. The arrangement of the fill is 
shown in Fig. 2.1. The crest zone was widened considerably 
after design conditions were revised during construction. 

The dam axis is curved upstream on a 10,000 foot 
radius arc, and located so as to optimise wedging action 
between the abutments. Under reservoir loading, this should 
give extra safety against tension cracking of the core. 

Surveying of the structure was based on a local 
coordinate system in terms of stations on and off-sets from 
the long chord. 

Adequate safety against normal and drawdown conditions, 
and also severe earthquake loads, was provided in the design. 


Basic features, some later revised, are shown in Table 2.1. 


Zoo HEStOory or Dam Construccion 


Twin 45 foot diameter concrete-lined tunnels were 
constructed in the left abutment between 1965 and 1967. 
The river was diverted in November 1967 after construction 
of the closure dykes, and cofferdam construction, riverbed 
excavation, core contact zone grouting and minor fill 
placement continued until March 1969, when the major fill 
placement program commenced. Dewatering of the construction 
area was by tube wells which were phased out when the fill 
height was beyond any danger of overtopping. 

Some foundation layers of fine sand in the riverbed 
exhibited non-dilatant behaviour under high stresses during 


testing. Thus, all fine sand was removed from highly stressed 
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Dam Axis 
Max. 2475 Ft eH gy 2500 F4 
Min. 2320 ft. \S— M3 
4 2 2 
TE. Sip Sal 
Upstream mas 
M2D Downstream 
C'dam 
ee sade. Orig. Riverbed ) y~ El. \885 is 
ete 
Blanket | Bedireck ay 
Grovtin 
200 ° 200 400. SECTION TYF 
Ss eee 
ZONE DESCRIPTION 
MI Core, glacial fill in 25cm (10%) layers. 
M2 Main shell, sand ana gravel in 30 cm (I2") layers, 


changed cbving construction to 45 cm (18”) layers. 


M2DI Inner zone of poorer M2 materials. 

M3 Core support Zone, sand and gravel or vock in 15 cm (6’) 
layers. 

M4 Outer shell, sand omd gravel or rock in 60cm (24”) 
layers. 

Drawdown 

Zone Gravel, cobbles and boulders or rock in 60 cm (24") 

layers. 

R.O. Original River Overburden 


FIGURE 2,1 FILL ARRANGEMENT AS DESIGNED 
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areas of the foundation because of the possibility of 
liquefaction under earthquake loads. Coarser overburden 
ana till was left in place except in the core contact area. 

Major fill placement occupied four construction 
seasons roughly from April to November each year, depending 
on climatic conditions, and was essentially completed by 
November 1972. The average placement properties of the 
materials shown in Fig. 2.1 are listed in Table 2.2, and it 
will be noted that significant changes were made after the 
1969 construction season. Design properties quoted for the 
dam materials are given, for comparison, in Table 2.3. The 
different fill properties have important consequences detailed 
later on in this thesis. A comment should be made concerning 
zone M2DI. In this zone the contractor was allowed to place 
fill which was available much closer to the site. These 
materials were generally inferior to zone M2. As a result 
of the instrumentation program it can be stated that in 
general the materials were of high quality and good performance, 
and given the precautions taken in design (Webster, 1970), 
that the structure is sound and successful. 

A concrete spillway structure, outlet works, and 
inlet facilities for future power development were constructed 
as fill was placed. An underground power station is at 
present under construction in the right abutment. An 
extensive survey network will monitor regional movements, and 


the University of British Columbia is managing a seismic 
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TABLET 2.2 


FILL PLACEMENT STATISTICS 1969 - 1971 SEASONS 


Year Material Yq (pct) wW(3) VY (PCr) 
1969 M1 fhe ES 8.8 143). 0 
M2 LAG, 1 Lao 158.0 
M2DI 144.4 7.6 156.0 
1970 Ml ates 6.3 150.4 
M2 146.6 fr | 1527,0 
M2DI 145.6 4.4 162i. 4 
1971 Ml eS: 6.2 T5150 
M2 146.5 pes) ES2n.0 
M2DI 145.6 4.7 1520 
1972 Not available 


1972 materials assumed to be comparable with 1971 
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array studyng the reservoir area. 
2.3 Design Revisions During Construction 


Placement Criteria. Following the first 
construction season important changes were made in the 
addition of fill. In the shell zones layer height was 
increased without significant change in compacted dry density, 
but the reduction in moisture content lowered the bulk 
density somewhat. Little change was caused to the material 
stiffness. In the core zones, the moisture content was 
reduced from near optimum (modified proctor) to 2% dry of 
optimum. Appreciable gain in compacted dry density was 
achieved and the core material became much stiffer than in 
the, 1969 season. 

The zone M2DI comprising "weaker" sand and gravel 
actually varied quite widely in properties, depending upon 
the source of the £111 at any particular level. Further 
discussion is made of this variation in Section 4.5 where 
the derivation of material properties is discussed in detail. 
For convenience this zone was treated as if it, too, was 


distinctly different after the 1969 season. 


Widening of the crest. When it was realized that a 
large slide or earthquake could cause a wave to overtop the 
dam careful consideration was made of the erosion hazard 
to the structure. Above elevation 2320 ft (707 m) the 


downstream slope was increased to 1 on 1.5, and the crest 
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width adjusted to 111 ft. (33.8 m), and details of the 
revised. crest geometry are shown on Fig. 2.2. With this 
arrangement erosion damage was not likely to cause general 


failure. 


2.4 Instrumentation Program at the Dam 

CASECO has produced detailed reports on the instru- 
mentation program for each year since the program started. 
Webster and Lowe (1971) provide a concise summary of the 
instrumentation of Mica Dam and comment generally on 
instruments which are available for monitoring high 
embankments. 

Several vertical or near vertical movement (MV) gauges 
provide the data on settlements, strains, and lateral 
movements which were extensively used in this investigation. 
Locations of these devices are shown on Fig. 2.3. They 
consist of telescoping sections of grooved casing installed 
with backfill in the embankment. Settlements were read using 
a field-developed latch-cone device or a radiosonde, while 
lateral movements were measured with a slope indicating 
instrument developed by SINCO of Seattle. Some details of 
installment of the MV gauges, and of the latch-cone device, 
are given in Appendix C. 

Despite some failures for various reasons, the 
piezometers functioned well. Earth pressure cells yielded 


disappointing results in general (Skermer, 1974). Strain gauge 
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28 
results were ignored until after construction machinery 
had finished working, due to "noise" produced by dynamic 
loading from the moving construction equipment. Settlements 
and incremental strain data from the MV gauges can be used 
with confidence, but due to the very small horizontal 
movements in the embankment during construction the horizontal 
movement data from the MV gauges must be interpreted in a 
qualitative manner only due to the limits of instrument 
accuracy. Despite some shortcomings, the program provides 
far more data than can readily be assimilated, and is quite 
successful and useful for monitoring the structure or checking 


analyses. 
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CHAPTER: *LiE1 
FINITE ELEMENT ANALYSIS OF EMBANKMENTS 


3.1 2D Versus 3D 

Prior to the development of electronic computation 
techniques, rigorous stress-strain analysis of continua was 
confined nearly always to a 2D simplification because of 
the enormous complexity of the mathematics. Even then, 
solutions were developed only for certain cases of over- 
simplified geometries. With the advent of the direct 
stiffness method, the finite element concept, practical 
matrix algebra and digital computers, great steps forward in 
analytical power were made. Even so, simplification to a 
2D state is still almost mandatory because of the logistic 
and economic problems of a 3D analysis. 

Clough (1960) first used the term "finite element 
method " in relation to 2D plane stress analysis. Using 
the 2D state very complex analyses have been made and in 
soils engineering developments have in many ways completely 
outstrippedy any) practical applications. Soils and rocks are 
complex materials, in general our limitations are not in 
analyses per se, but in understanding and being able to 
measure basic behavioural characteristics (if these are 


recognized) of the materials. Finite element 2D analyses 
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30 
are now capable of answering most of the questions that may 
be placed before them, and a matter of increasing concern 
is) whether the 2D’ result is good enough in a 3D world. 
Experience suggests it is, in most cases, but only by 
evaluation of the two approaches can the question be truly 
answered. 

2D programs can be run on modern computer facilities 
at a very nominal cost, considering their power and impli- 
cations. With 3D programs, a very large sophisticated system 
is necessary to run even a modest analysis at a cost which 
may not be justifiable in terms of the additional information 
obtained. 

The isoparametric finite element concept, powerful 
and generalized in 2D or 3D, has been widely adopted because 
of its versatility (Zienkiewicz, et al, 1969). Clough (1969) 
evaluated the most suitable element forms for 3D analysis. 
The most advanced system available at the time of writing, 
to the writer's knowledge, is documented by Bathe, Wilson, 
and Peterson (1973). In geotechnical engineering the most 
practical 3D results reported were provided by a program 
developed by Krishnayya (1973a). 

Recent studies (Palmerton and Lefebvre, 1972; 
Eisenstein, Krishnayya and Morgenstern, 1972a) have clearly 
indicated conditions in which 2D analyses do not agree well 
with 3D analyses. The 3D situations are primarily a dam in 


a narrow, steep-sloped valley and dams with a thin central 
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a 
or inclined core of markedly different properties from the 
shell. Mica Dam satisfies both of these. The increasing 
development of high fill dams in narrow and irregular 
valleys accentuates the need to know the 3D effects in finite 
element analyses, in order to more adequately predict 
ammportant *features*®suchvasicrackingsor~hydraulicefracture. 

Previous work (Covarrubias, 1969; Strohm and 
Johnson, 1971, Eisenstein, Krishnayya, and Morgenstern, 
1972 a,b) has provided good broad outlines as to factors 
intr luenci ng?’ cracking .oThercontrol oficrackingninéeearth 
dams is of paramount importance and is possibly the prime 
use for finite element analyses which can successfully predict 
movements and stresses. It remains for a predictive 3D 
analysis of a real structure to be reported, and this thesis 


work may assist such developments. 


ce Representation of Soil Stress-Strain Behaviour 


The classic formulation of stress components in terms 
of strain components for an ideal, isotropic, elastic solid 
can easily be written in terms of the familiar Hooke's Law, 
using Young's modulus and Poisson's ratio. Young's modulus 
is a measure of normal-stress/strain response, while Poisson's 
ratio is a measure of lateral deformations associated with a 
particular directional deformation. Neither of these 
traditional parameters has any clear significance when we are 


dealing with stress states in a continuum. The Bulk and 
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87 
Shear moduli, which can be derived from Young's modulus 
and Poisson's ratio (or vice versa), do have a physically 
separable significance »with,regard ito,octahedral normal and 
shear stress states respectively (Domashcuk and Wade, 1969). 

Natural materials exhibit non-linear, anisotropic, 
stress-or-strain-path dependent stress-strain response. 

Only by the crudest approximations can they be treated as 
ideal, isotropic elastic materials, and often the stronger 
Materials can be realistically represented in this simplest 
way under normal working stress conditions. Finite 

element techniques relate stress and strain fields through 
constitutive equations, usually based on the simplest forms 
of elasticity. These latter result in easy fast calculations 
which can yield good results in general terms (e.g. 
Covarrubias, 1969) but which have limited success as flexible 
predictive tools. The other extreme is the 21 independent 
coefficients relating 6 stress and 6 strain components at an 
infinitesmal point, all of which may vary in a complex way 
with material and with stress state. 

Volume change due to shear stress alone, or rotation 
of principal stress axes, are one common deviation from ideal 
elastic behaviour in soils. Another is the incorporation of 
plastic deformations into elastic response. Unhappily, the 
present state of soil testing is such that only very limited 
data for an analysis invoking complex behaviour are available. 


Smith and Kay (1971) describe a more satisfying approach to 
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39 
soil deformation behaviour, but one which is still beset 
with problems of inertia of the discipline and mystifying 
appearance to most of the potential users. 

The most successful approaches yet made to soil non- 
linearity assume piecewise linear incremental elasticity. 
Footings have been analysed using a Young's modulus derived 
from Triaxial tests, and an assumed Poisson's ratio 
(Radhakrishnan and Reese, 1969). Stresses and movements in 
Oroville Dam have been calculated using iterative techniques 
based upon Kondner's (1963) proposals for hyperbolic stress- 
strain relations for soils (Kulhawy and Duncan, 1970). Data 
for non-linear variations in Young's modulus and Poisson's 
ratio are obtained from triaxial tests and approximated by 
hyperbolic functions. kKrishnayya (1973a) analysed cracking 
in Duncan Dam using Triaxial test data stored in stress 
invariant form, but limited Poisson's ratio to a maximum 
of 0.49. Penman and Charles (1973) assumed values of 
Poisson's ratio and used the concept of constant equivalent 
compressibility from 1D behaviour to account for changes 
in modulus with stress level. 

As pointed out so rightly by Smith and Kay, the 
results of such analyses are surprisingly good considering 
their deviations from the real stress-strain response of the 
soils. It may well be that simple methods give really good 
results for typical embankments, and until time effects or 


unusual stress paths or failure is considered, there is no 
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34 
need for further complexification, just a hard-headed sense 
of the realities and limitations involved. 

In this investigation of Mica Dam the Bulk and 
Shear moduli formulation, as proposed by Clough and Woodward 
(1967), is used. The effect of modulus change with stress 
level is introduced by allowing the One-dimensional modulus 
to vary with octahedral normal stress, and a generally 
constant value of Poisson's ratio is assumed in calculating 
tables of K and G versus o for a material. This method 


Cir 


is more fully described in Section 4.5. 


3.3 Shear Failure Behaviour of Soil 

Kulhawy et al. (1969) suggested that after shear 
failure of an element further stress-strain response could 
be modelled by keeping K at its pre-failure value and 
decreasing G towards zero. kKrishnayya (1973a) computed shear 
failure using the same criteria as Kulhawy et al, and reduced 


the Shear modulus to 


ae 


50 after failure. 


G = 


The value of K was determined by unconfined compression test 
data, when he used the stress-invariant triaxial data approach. 
Otherwise, as suggested by Kulhawy et al, K keeps its pre- 
failure value, and this is the approach used during the 
present investigation. Stress redistribution and plastic 


deformation are only partially represented by this formulation 
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(Smith and Kay, 1971). Consequently the results of shear 
failure in the reported analyses of Mica Dam should be 


viewed with caution. 


3.4 Pore Pressure Effects 

During fill construction pore pressures are generated 
by total stress changes and dissipated through drainage. 
Particularly in thin cores of silty clays orgtacial tilis 
there may be appreciable dissipation over the construction 


period. Koppula (1970) utilized the governing equation 
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2D coefficient of consolidation 


where a 


u =repore water pressure 
Y = bulk™density or "the i111 
= = Au (Bishop, 1954) 
B aon 
Ao 

i 
oh -= time rate of increase of fill height. 
at 


which can be applied to vertical cores quite accurately. 
Krishnayya (1973c) developed a 2D finite element program 
which solves this equation, and which will be used in 
conjunction with total and effective stress analyses in 
another part of the Research Project involving Mica Dam 
(Law, 1974). Chan (1971), applied the equation in finite 


difference form to simulate pore pressures in Mica Dam, 
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quite successfully except for the effects when foundation 
tube wells were deactivated, which he could not allow for. 

the dafficultyin stress analysis is dan simulating 
the effects of drainage, as the deformation parameters can 
vary quite significantly. Partial drainage can be dealt 
with by doing analyses for the two limits of no drainage 
and full drainage (e.g. Chang and Duncan, 1970). The correct 
stress-strain relations can be determined from laboratory 
tests that simulate the proper stress paths and partial 
drainage consistent with field conditions. This can be 
impractical. In general the most successful attempts use 
consolidated-undrained test results as a reasonable medium. 
Fully drained conditions correspond to consolidated-drained 
tests and undrained conditions to unconsolidated-undrained 
tests. Poisson's’ratio will approach 0.5, the condition 
of incompressibility, only for total stress analyses under 
undrained conditions. The importance of the best possible 
simulation of drainage conditions has been emphasized by 
Duncan (1972) and Lowe (1972). 

In the present work, the fill materials were all 
considered as free draining except in the core. Thus 
there was no distinction between effective and total stresses 
in these materials. In the core, laboratory tests indicated 
that only in the 1969 material were significant pore pressures 
likely to be generated, and a reasonable judgement was made 


to obtain representative moduli (Refer to Section 4.5). 
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a) 
3.2 Pneremental Anal ysiSVof Fill *Construction 

Clough and Woodward (1967) demonstrated that 
incremental finite element analyses are required to calculate 
displacements in an earth fill that would be measured in the 
field. For a homogeneous dam, they indicate that stress 
distributions from a simpler "gravity turn-on" analysis are 
reasonably close to those of the incremental analysis. Since 
actual ‘constructton"involves” Layers*fiar too thin teebe 
individually applied in a finite element analysis, investi- 
gations have been made of the detail required to gain 
respectably good results with the finite element approximation. 

Kulhawy et al (1969) discuss the layering problem 
in considerable detail, and suggest a means for determining 
the minimum layering required for a suitable accuracy. Other 
factors such as overall mesh size, computer resources, job 
requirements etc. can be expected to play a role in economic 
decision-making. Naylor and Jones (1973) used a remarkably 
coarse mesh in a separate analysis of Llyn Brianne Dam and 
construction was simulated using only 4 layers and a high- 
order element. They report results as successful as Penman 
and Charles' (1973) detailed analysis. 

Other, more subtle effects arise with layering 
analysis. To what extent should shear between placed material 
and new fill be considered? Clough and Woodward used a 
"reduced modulus" to allow for reduced shear interaction 


between a new layer and existing fill. However, several 
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techniques ave usedsiniconstruction toxkey aenewflayex: into 
the existing surface, particularly with more plastic materials 
which can actually shear to residual under construction equipment 
loads (Sherard et al, 1963). Another approach is to assign 
the properties of a dense liquid to the new increment of 
fill in an analysis. This is neither necessary nor 
realistic. 

Two fundamental approaches seem to have developed in 
incremental analysis. It is useful to examine them for their 
own sake, although it is doubtful whether their results 


would be significantly different. 


Approach (A): 


1. The first layer of elements is placed and assigned 
stresses so that OF is the overburden pressure 


at the centroid of an element and 05 and 0, are 
calculated using the coefficient of lateral 
pressure (which may be calculated from Poisson's 
Hat io). 

2. The structure stiffness is calculated and nodal 
loads are applied to the upper surface corresponding 
to the weight of the next layers ovIf mecessary 
these steps are repeated using the "average stress" 
approach. 


3. The next layer of elements is added and stresses 


assigned in these as before. The stresses and 
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displacements due to the next layer of fill are 
then added to the previous results. The process 
iswrepeated until all fill loads have been 
added. 

This approach is exemplified by Kulhawy and Duncan (1970) 


and Palmerton and Lefebvre (1972). 


Approach B: 


1. The first layer of elements is added and stresses 
and displacements due to self-weight of this layer 
are calculated using the layer stiffness. 

2. The next layer of elements is added and stresses 
and displacements for this structure are 
accumulated with previous results. If necessary 
the "average stress" iterative approach is used. 

3. The process is repeated until all fill has been 
placed. Thus nodes along the top of the fill have 
displacements due to the self-weight deformation 
of the last-placed layer. This situation does 
nobiloceur ioeiplacing <and!imonttoring of field 
measurement devices. Hence, the construction 
displacements of nodes on the upper surface of 
the nth layer are computed by subtracting from 
their total displacements the displacements due 
to self-weight of layer n. 


This approach is exemplified by Krishnayya (1973b) and is 
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the method used in analysing Mica Dam using the program 
FENA 3D. Intuitively, for an infinite number of layers, 
this process would be the more correct. In the normal 
Situation it is unlikely that there will be significant 


difference between the approaches. 


Theoretical Studies Using Program FENA 3D: Eisenstein, 
Krishnayya, and Morgenstern (1972a) compared results of 3D 
analyses and 2D plane strain analyses of the longitudinal 
section of a dam ina 1:1 valley. The results are 
summarized again here because of the light they shed on 
results of the 3D Mica Dam analysis. All the analyses used 
linear properties and varying numbers of lifts. 

1. There should be about 10 lifts to compute strains, 
stresses and displacements to agreeable accuracy using a 2D 
analysis. For a 3D analysis economic considerations may 
reduce the number of lifts, but this reduction should be as 
small as possible. Refer to Fig. 3.1. 

2. Only for homogeneous embankments will a plane 
Strain analysis of a longitudinal section in a steep valley 
yield satisfactory information about stresses and displacements. 
The effects of different stiffness in core and shell force 
us to use a 3D analysis to obtain satisfactory information 
(Pig.eg.2). letebvre et al. (1973) discuss in more detail 


the effects of 2D and 3D variations in homogeneous embankments. 


veseqal ep ae peg way ae Ai 

Leayor ait m1 Lsoerius 108 aj.edbLuow seezowg ait 

gnepitiapte ed [liw otend send plotiiow ei cf aolseusie 
acdsnihiie-aaa-idenied mma 

sony 6FvRe Sa 


Wiegeneata 6 3de _IMg_mstgo18 potet ealbus2: tageggeeess 
GE 70 adiuaer ABIAgHOD AARTELT ameawneeteM Bas CREE, 
rsniSutipdol sit Yo aseyféns nisiew snsiq Os bas sonyens 
ats esluaar oat Lettay: Gh « ot oabim 20 molaomg 
io heta yohs sdoki og to sausped sia ehepe sf 
boay esevinns ote ILA -.ciayions mG soiM GE eld to esitand 
_edlil lo ezedmun podyusr Bb eettaeqoug tssakl 
\2hisige. Stugwiop) op 235551 Of tuods sd blyeie eusd® .2 a 
di s beteu yosTUo0s eidsestps of etvemoniqekl bag seesezte 
¥en ahoitsisbianos oimonocs eieylens Cf 6 x07 valeyisns 
as ed Bluode solsouuexr eins Jud .224if to sadmun ond soubst 
1.£ 2obT of ishe® sefdigeog se Liem - 
Sreiq = ILim ainsiblasdm 2vcenayoust sokgia® -& =) 
follev quase se al. notisse Jans6utipool s te aseyions atsitse 
.eonememsigeib bas agzesxte tuode onl tesrotnd grodeetelsse bioly 
enyet Ife tvié sy02 ai ensattidz ingrettub 2o atdette oft 
AOk¥emoINi. ysoddRtet tea) abstiio ot ecegians mais eeu of Bo 
Liaiad axon mk ago2eib (EVEL) ss 26 ceeibasd RED 


Pane 3 
SOIL: E= 200KG/CM 
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FIGURE 3.1 EFFECT OF NUMBER OF LIFTS ON SOME IMPORTANT 


PARAMETERS (after Krishnayya, 197 3a) 
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FIGURE 3.2 COMPARISON OF HORIZONTAL STRESSES ALONG A DAM CREST 
FOR 2D AND 3D ANALYSES (after Krishnayya, 197 3a) 


CHAPTER IV 
FINITE ELEMENT MODELLING OF MICA DAM 


4zinDiscussivons and Site’ Visit 

Discussions at various stages of the work described 
in this thesis always involved some aspect of the dam's 
character. This was felt to be invaluable to the writer 
as, over a period of time, this constant reference to the 
structure resulted in a confidence and familiarity with data 
which it is felt is very important when the time comes to 
discuss results. Comment on this aspect is included here 
because the writer regards it as a very important phase of 
finite element analysis in general. Particular emphasis 
must be placed on discussions with CASECO staff. A visit 
to the damsite was made in September, 1973 where features 


of the project were seen and discussed at first hand. 


4.2 Construction and Instrumentation Drawings and 


Report 
CASECO Consultants provided the following information 


on the dam structure and materials: 

(a) Field Instrumentation reports for the 1969, 1970, 
and 1971 construction seasons, and a set of preliminary 
drawings from the 1972 report. 

(b) Reports of triaxial tests performed by other 


agencies; of the core, shell, and river overburden materials. 
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(c) A selection of construction drawings showing 
the planned and actual excavations which occurred in the 
core zone. 

(dj) A series.of fill progress drawings showing the 
transverse profiles across the valley. 

(e) Sundry details including some finite element 
analysis results, and construction drawings for the upstream 
and downstream cofferdams. 

These resources were continually employed in 


determination of structure geometry and material properties. 


4.3 Synthesis of Mica Dam Geometry 


A master plan was prepared at a scale of 1 inch = 
100 feet on which was reproduced the detail of all available 
construction information, the result being a contoured 
representation of the site after all overburden removal 
and core zone excavation was complete, just prior to the 
placement of fill. 

Since the prime concern in developing a model of 
dam geometry was to successfully represent the slopes and 
slope changes of the abutments, a profile was prepared (from 
the master plan) of the intersection of the longitudinal 
central surface and the valley surface. (These terms are 
explained in Fig. 4.1) From this profile the important changes 
of slope were identified, and transverse abutment-surface 
profiles constructed at these points. Fill profiles at each 


of these locations were then added, so that the dam geometry 
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PLAN VIEW 


(Outlines after Skermer, 1974) A TYPICAL TRANSVERSE 
SECTION B-B , NORMAL 
TO LONG CHORD. 


Va lley 


Surface 


Valley Surface 
Trace of Longitudinal 


Central Surface at this 
Transverse Section 


The surface swept out by the sectioning lines A-A, in moving 
from station to station along the long chord, is termed 


the LONGITUDINAL CENTRAL SURFACE. 


FIGURE 4.1 EXPLANATION OF VALLEY SURFACE AND 
LONGITUDINAL CENTRAL SURFACE 
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46 
was described by a series of vertical sections transverse 
to themtonge chord. 

The element mesh could then be assembled by sub- 
division within the vertical sections, having regard to 
abutment surfaces and internal structures Of@tiecwriime At 
this stage the two-dimensional visualization process broke 
down, and a model was constructed. The series of vertical 
sections were reproduced with coloured marking pencils on 
sheets of lucite, which were mounted in their correct positions 
on a base-board covered with a copy of the master plan. This 
arrangement, giving visual access to the 2D and 3D geometry, 
was invaluable at this stage and even more valuable at other 
stages, and the writer emphasizes that such a 3D model is 
absolutely essential for an analysis of this sort. Even 
so, a great amount of effort is necessary to mentally visualize 


small details. A view of the model is shown in Fig. 4.2. 


4.4 Subdivision Into Finite Elements 

The inaccuracies of poorly shaped and skewed elements 
become crucial at this point, as too large a subdivision will 
be too crude, too small a subdivision will be impossibly 
expensive, and anything between these limits must have 
satisfactory elements and representative geometry. The 
process requires considerable experience and since the time 
involved is measured in weeks or months, the writer emphasises 
that herein lies a cost at least as significant as computer 


costs. 
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A mesh regarded as "satisfactory" was finally 
developed. After considering the shape of elements there 
resulted a finite element mesh of 254 elements and 276 
nodes, yielding a structural stiffness matrix of 828 
equations with a semi-bandwidth of 183. After removing the 
boundary nodes (where some displacement condition was 
specified), there were 438 effective degrees of freedom. 
This mesh consisted of five lifts: 

1. river overburden after excavation 

ZeectOeo. gti) placed an the 19169 to 1972 
construction seasons respectively. 

Since one element height thus represents one full 
season of fill placement this is the coarsest mesh possible. 
Even thus, this mesh is as large as could be conveniently 
used for repeated calculations, and so is the finest mesh 
possible. 

The complexities involved meant that simulation of 
river overburden excavation was not attempted. Since these 
materials are probably highly overconsolidated this is a 
deficiency of the analysis. The results obtained for river 
overburden material were thus treated with extreme reserve, 
and in this area the worst element-geometry problems 
compounded the doubts. Further discussion on this topic can 
be found in Chapter V. The problem of the river overburden 
elements may have had significant effects on all calculations 


made, but it was felt that even so the overall results of 
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analyses using this mesh would be quite realistic and 
representative. 

During the process of subdivision several nodes 
were located between sections, using model and master 
plan. Extra elements were included to assist in modelling 
the river overburden materials. The mesh geometry is 
presented in Appendix B. Figure 4.3 shows the section at 
station 22+50, the main instrumentation plane and the 


assumed plane strain section. 


4.5 Derivation of Material Properties 

The program FENA. 3D contains an option to use 
stress-strain data from triaxial tests to interpolate elastic 
parameters for an element based on its stress state. This 
is an elegant but complicated and expensive way of performing 
the analyses. Accordingly most analyses were carried out 
using “engineering judgement" to select the "representative" 
linear elastic moduli. As pointed out by Skermer (1974) it is 
very difficult to obtain laboratory data which match field 
behaviour for soils, because of compaction and grading 
problems chiefly. In most cases incremental stress-strain 
data from the MV gauges were used. Laboratory tests by T. C. 
Law which are part of further research being undertaken at 
the University of Alberta, were available for core material. 
Reports of triaxial testing for CASECO were also available 


for some materials. A summary of available stress-strain 


: eset un 

7M os o 
debi xeving, gia ivi bd sete: 
‘ =- iene Sees 
tstanm bis tabon a : 
gikiGbow A arte, oF snk aa 
ai ssbu0miQ dae ath olptxedinm Gulandave 4 
th moLjasp. ons Aworle €,% sxvpit A xi baeqga 
az Das aiéia nol ssacemustand kom ad? : 


_ 


* 
— 
i 
oa 
Le 
tx 


ae 


® 


_ 
: 


4 


- 
; 
- 
> 


a 


sijesis seBfoqissn! of ajser Laixnias ms ht 
aifT .sdete ee9-t2, asa on -beged sHOmeS oe x0 paren 


pibmoarisg i6.Vvew oVLeqequs Sas 597657 Leniea ud 3s 


ai 32 (b90L) semmed? yd tuo besriog BA _tLubom 5 
blah? odem dostin mab yuothzednt nintdo af : 
enthery brs aeisnaqmos to seuplsd aoa 
nietie-sasxje [siremezont esses deba br Rtas z 
aed Need aseeet eaigisodat bean scow eppwtin’ UN afta 
$8 aedwenbny ind. waves: wektx03 ve sg on 


50 


OOL! 


0o6! 


oold 


o0EZd 


0oSZ 


(14) AatTa 


HSAW LNAWH TH ALINIA 
°06+7Z°VLS “NOILOES ASUHASNVUL €°7 HAUNOIA 


(14) dYOHD DNOT Wesa BNVISIQ ASUSASNVYL 
090! + ) 


9oo!- 


NOQOLVI01 ADNWD AW — 


(o13) €°0N 1411-€) 


AadON @ 


(72) DAQH>D abs MONA So4ATaIC APSVEMART 


‘O@4CS ATe2 ,MOLTOSe ZRAaVeMART Eve BAWOTA an 


_ : = Noam SWELLS ATLMEY - 7 


51 

TnHroOrmMatron eB Sigiven im Table! 4a. 

An approach towards non-linear behaviour of the 
soils was made by allowing moduli to vary with stress level 
in the elements, as described hereunder. As previously 
mentioned, elastic, isotropic materials possess two independent 
stress-strain parameters. The usual form is Young's 
modulus (E) and Poisson's ratio (uw). An alternative 
formulation uses the Bulk modulus (K) and Shear modulus (G). 
In this work the One-dimensional modulus (D) is ubiquitous. 
In general, all these parameters vary in their own way with 
variations in the stress state. Assuming isotropy and linear 


elasticity the following interrelationships exist between the 


variables: 
LD) (ee eae 
pee or aos 
E 
= "i 2(1 + w) 
2 
mee | DUR ue) 
ea) 
Poissons' ratio: In general Poisson's Ratio will vary 


depending upon the stress state of the soil. Computation of 
horizontal movements and minor principal stresses is more 
sensitive to variations of this variable than settlements and 


major principal stress (Skermer, 1973; Kulhawy et al, 1969; 


Covarrubias, 1969). 
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Skermer (1973) indicates the change of Poisson's 
ratio with stress state during an analysis, where non-linear 
triaxial data was used, and the principle is generally 
accepted. Results from the Mica 3D analyses (linear 
properties) indicate that the computed stress state for a 
Material does not necessarily agree with the relation developed 
between oF and Ot in Section 1.3. What deviation was due to 
stress transfer or computer roundoff errors, and what deviation 
may be an indication that the material should have exhibited 
a change in value of Poisson's ratio could not be determined. 
The deviation shown in Fig. 4.4 is for shell zone M2 material. 
The writer argues that the assumption of a reasonable value 
for Poisson's ratio is a satisfactory first-order approximation 
to the true behaviour. 

One-dimensional compression: Particularly at the 
early stages of construction any lateral movements of a fill 
are negligible compared to settlements. It is satisfactory to 


assume one-dimensional compression for which the incremental 


modulus 


Povere 


Aevert 


This one-dimensional modulus, as shown earlier in this 
section, is nothing more than one of a set of two required 
constants for simple linear elastic analysis. It perhaps 
best describes the stress-path undertaken by a soil element 


during fill construction, particularly of high sais lor Stitt 
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mMaterials@such®asi MicalDamelradnithi s structure) lateral 
movements were very small and field measurements of the D 
modulus satisfactorily reflect the stress-strain behaviour 
during construction when used with an appropriate value of 
Poisson's Ratio. This modulus was invariably higher than 
equivalent moduli calculated from triaxial tests on the fill 


materials 


Allocation of material prOpereies= (Onwene basis 


of field and laboratory data six different stress-strain 
modes were selected for the Dam materials. A seventh material 
was allocated due to the higher bulk density of the M2 shell 
material in 1969. The materials are: 

1. River overburden 

2.8 onell M2) (sand; and Gravel) 1970-1971 

Speemcore Ml (Glacial bulb) Lo 70-1972 

4. Inner Shell M2DI (Sand and Gravel) 1970-1972 

5. Inner Shell M2DI (Sand and Gravel) 1969 

Beaeee ore: Mil (Glacial. P21i).eL959 

7. Shell M2 (Sand and Gravel) 1969. 

Linear elastic parameters were selected as described 
hereunder. In all cases u was chosen as a reasonable value, 


the most important influence being the relationships: 


K = le Sino’ (Brooker and Ireland, 1965) 


K 


and U = iE (linear elastic theory) 
fe) 
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Nonlinearity was introduced by allowing the moduli 
to Varyridepending upon stress level in the elements. - It 
was assumed that the One-dimensional modulus varied with the 
octahedral normal stress in the element. From the field 
curve of vertical stress versus vertical strain a plot of 
the D modulus versus Oe was made. Using the relationships 


ert 


of Soa pes and OF and O71 and OP et 


material by averaging the results of linear analyses of the 


obtained for the particular 


dam , D was expressed in terms of om Assuming appropriate 


Cree 
(usually constant) values of uw a table of K and G moduli 
for various Crsestee levels was assembled for computing. As it 
turned out the theoretical relationship discussed earlier 


between o, and che was sufficiently good to represent the 


1 qe 
finite element behaviour, so that a table of K and G versus 


o, would have been just as appropriate. The process of 


a3 

obtaining the multilinear properties is shown in Fig. 4.5. 
The multilinear approach thus enables the stiffness 

of an element to vary with stress state. It does not take 

into account dilatency of the materials nor the variation 

of Poisson's ratio in general. The stiffness is derived 

from one-dimensional stress-strain data and has similarities 

to the approach described by Penman, Burland and Charles 

(1971) and Penman and Charles (1973). In the program FENA. 3D 

an "average stress" process (as described by Kulhawy et al, 


1969, and used by Krishnayya, 1973b) was used to obtain the 


best stiffness at each layer of analysis. 
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Gauge Data Curve: 


oe 


€ 
OD etBortoe! vert 


plow. 


Results of Linear Analysis: 


(1) Se rent 


at gauge location. 


IDs relationship 


(2) usually ¢€ € 


vert 1 


Convert Gauge Data curve 


to oO /e, curve. 


Compute Slopes of o7/e, curve 


at various values of 7: 


Plot D modulus /O curve. 


Results of Linear Analysis: 


eee relationship with et 


c 
at representative points. 
(2) Average to obtain Mer !P7 


curve, 


Derive D/O” curve 
oct 


Oedometer Test Data: 


oO ./e curve. 
pee 


Assume suitable values of MM 


at©O” | values on D/C 
oct oct 


curve. 


Calculate K/C~ and G/O~ 
oct OCE 


curves and represent as a series 
of linear segments for computer 


tabulation. 


FIGURE 4.5 MULTILINEAR PROPERTIES SELECTION PROCESS 
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Properties used in the linear analyses are shown in 
Table 4.2 and the multilinear properties are shown in Table 
Ai. 3. 

Core material: By coupling the results of high-pressure 
oedometer tests and all-round pressure tests it is possible 
to obtain two independent stress-strain relationships and 
hence to calculate two independent elastic parameters which can 
vary with stress state. Such work is currently being 
performed on the Mica Dam core material as part of further 
analyses within the Research Project, by T. C. Law (1974). 
While it was possible to incorporate this extra information 
into the 3D analyses, a decision was made to use simple 
judgement in assigning an essentially constant value of uL. 
This simplifies the selection of data for computation but 
remains to be justified by comparison of analytical and field 
data. The 3D analysis was felt to be too coarse to justify 
any more than simple approximations to the truly complex 
nonlinearity of soils. Nonlinearity is most evident in the 
finer grained core material, so in this region the greatest 
simplifications have been made. Some oedometer test data 
are shown in Figs. 4.6 and 4.7, and this one-dimensional data 
was used to select D moduli. 

Shell zone M2: Field data was available from several 
MV gauges but MV11 and MV12 were deemed to give the best 
results, being influenced essentially by valley arching only. 


A plot of data used to calculate D moduli is shown in Fig. 4.8. 
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Shell inner zone M2DI: CASECO were able to supply a 
Study ofaehe Variation of materials in. this sone. who 
as described earlier was comprised of poorer materials from 
different sources. Their plot of the variationsoiamodud 
and materials is shown in Fig. 4.9. Skermer (1974) was able 
to model the distribution of these materials with greater 


Certain his shiner 2D. mesh. 


River overburden: Triaxial tests were performed on 


this material. Field settlement data indicate a bulge at 
the overburden - shell interface, suggesting some disturbance 
during construction works. In general it is to be expected 


that these materials would behave as very dense, over- 
consolidated granular fill, and their insitu deformation 
properties would be correspondingly much stiffer than measured 
in a normal test. This factor was taken into account when 
selecting moduli, as the field data is necessarily quite 
limited for the river overburden. A plot of some gauge 


dita 25 Siiown in Fig. 4.10. 


4.6 Synthesis of Bedrock Movements 


A tremendous volume of river overburden material 
Was excavated for the core trench of sthe dam. » Field 
measurements indicated that there was an unexpectedly large 
heave of the bedrock at this time, and that some heave 
actually continued as fill was placed during 1969. Rock 


settlements in this zone were greater than expected, probably 
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due to loss of tightness in the rock mass due to heaving. 
Law (1974) correlated some bedrock heave with pore pressures 
in the rock, but not enough data is available from the 
excavation period to make conclusive arguments about the rock 
movements. 

Krishnayya (1973a,b) analysed cracking at Duncan 

Dam, which was due to excessive differential foundation 
settlements, by adding observed foundation movements at 
appropriate boundary nodes of the finite element mesh. It 
was very easy to adapt the data on bedrock settlements at 
Mica Dam. By plotting the settlements at pertinent sections 
of the dam the incremental settlements for each lift of the 
analysis were determined and applied at the boundary nodes. 
Some interpolation and "guestimation" was required but in 
general the bedrock settlements are relatively very 
small and a representative pattern is a fairly simple process. 
For the larger settlements there was some concern as to 
accompanying horizontal rock movements and after some very 
simple analyses of rock slopes under fill load a rough guide 
was obtained, depending upon rock slope and position under 
fill. The horizontal movements were most likely too small 
to be influential and in retrospect it would have been simpler 
to ignore them completely. longitudinal and transverse sections 
showing bedrock movements as measured at MV gauges are shown 


in Figs. 4/11 and 4.12 respectively. 
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TABLE 4.3 


MATERIAL PROPERTIES FOR MULTILINEAR ANALYSIS 


Octahedral 
Normal Stress u D K G 
Level (Ksf) (KsE£) (KsE£) (Ksf£ ) 


River Overburden (1) 


Ore0 On 7800 4950 2154 
0.80 oye ik 8500 5385 2348 

steve 1, Spe ll 9800 6210 2708 
1S egy ieee le LEUCG 7040 3066 
36.80 oS is 12300 7790 Gye js 
48.40 | Cou 13000 8250 3590 
very large ae gil 13600 8620 sT55 


Shell Zone M2 (1) and (7) 


0.0 G25 10140 6020 3100 
4.50 Cae 10000 530 3055 
Ney Pe Le) eee 9050 5360 2766 
41.00 Vezo 9500 5630 2900 
very large 0.28 10000 5230 3055 


Core Zone Ml 70-72 (3) 


De Au) Usoo 1695 dd Wea 431 


egal) ae 2220 1456 548 
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TABLE 4.3 (Continued) 


Octahedral 
Normal Stress u D K G 
Level (Ks£) (Ks£) (Ks£) (Ksf) 
FLO 0.38 2200 1456 548 
SE RSW 0. 38 6000 3970 1523 
44.50 0238 8160 5400 2070 
very large 0.38 11000 7200 2792 


Shell Zone M2DI 70-71 (4) 

0.0 C222 7350 4480 2186 
very large 0722 7350 4480 2186 
Shell Zone M2DI 69 (5) 

0.0 0.29 8400 5120 2500 
very large 0.29 8400 S120 2500 


Core Zone Ml 69 (6) 


g.0 0.32 214.3 G0. 0 Aled 

i 67 0. 37 3.0.40 2267.0 64.2 

a. 70 0.35 6.0 480.0 146.3 
21657 0 0.36 i272 905 278 
19.30 0. 35 2000 1390 462 
29550 Vass 3500 2427 808 
42.50 O@35 5200 3600 1200 
very large 0.35 5200 3600 1200 
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TEST DOT-2°: standard proctor compaction, gradation |, -%" sizes 
w= 8.95% ¥4 = 132.1. pf. 
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linear analysis. 
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OF several tests, this one most closely followed field measurements 


FIGURE 4.6 CORE ZONE M1, 1969 MATERIAL: 
OEDOMETER TEST (after Law, 1974) 
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TEST: DpoT-6 kneading compadion, grad. | -%4 sizes, PLC MEMES leek 


DOT-14- modified proctor ; S24 7 ee 


O, | total and eftective 
(Ksf) principal Stress 
ie (negligible pore pressure 


generation ) 
| 
averaged 
100 modulus 
Ichosen for | 
| linear analysis 
90 | 
L other + 
Be | not as DEORE vy 


| of field conditions 


70 / 
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40. Working 
Stvess 
range. 
30 
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(0) 
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oe 
principal strain cS (7%) 
FIGURE 4.7 CORE ZONE M1, 1970-72 MATERIAL; 
OEDOMETER TESTS (after Law, 1974) 
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MV _ Gauge Data for Shell Zone M2: 


stress transfer effects probably minimal, 
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FIGURE 4.8 SHELL ZONE M2: FIELD MV GAUGE DATA 
(after Law, 1974) 


66 


(OONSVD) (OO+EZ°VLS) NOILINAIULSIG TVIYALVW OL GauVvdWOd 
(OS+7Z°VLS) OLTAN WOU SHILYAdOUd 


67 


>TG@W ANOZ WANNI TISHS 6°7 ZaNOIA 


(eee 


‘NadeNasano BaAizt 
3NI1 NONVAYOx3 


SJIN3BNTSANOD | oosi 


| 
= = j 
eee acl —f~_pabruanp +9A5 ulnsys jor1paA VANS doom Slr % 
be Si Nd JOP 1A 
JO} snjnpow abovann IIe Ea aa Bi 33IH) ATLL +8 


33D TILL «62 
Bary doom 1 -y| 006! 
Zanaasnod ec 


AINSBWdIWYINE SC | 


yanrd doom ov], 
gam atin szt% 
SONSMANOD OF 


‘gyva Jafio, wos, ‘aah yooe2 
uy _g sAjnpow abosany 


*SINaNMTANOD 
Guy ‘LNaWdWYONa oolz 
‘gana qoom AINIvVUW 


OLer 
‘Suosvas {1-016 


Jo} Q SA|Ppow abosano payybiam 4aiW> SALLI 


KINIWW (4) 
‘}Aan F 
aa SD VayOL ‘q NAS 


(354) G cooe! coos! 000Z} 0001) 00001 0006 0008 cool 0009 ANOZ 3O dOL 


a——_= 
en 
- 


a ee er os Ss eee 
e 1 
Se eae na?) SN ae ha re ee 


) taw> te th 8) Sr We et. ee «Bo ee 


(Oc+4.5 ATS) OLVH WOSTF SITIAA(OFF +: 10SM BMok SPW Jae oe. 4 Bee ae 
(OG8CAD) (OO4E5 ATS) AOITUAIATEIO IATASTAM OT G3SAaMOD 


120 


ners 


110 
(Ksf) 


100 


80 


70 


60 


20. 


68 
MV Gauge Data for River Overburden: 


styess transfer effects to be alloweol for. 


best f+ curve a better choice for 
to field data from deat Me 
MV gauges moclulus. 


average modulus 
Chosen for 
linear analysis 


ot data poivits from MV Qauges 


in River Over burden. 
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FIGURE 4,10 RIVER OVERBURDEN: FIELD MV GAUGE DATA 
(after Law, 1974) 
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CHAPTER V 


COMPUTATIONS USING THE THREE-DIMENSIONAL 


FINITE ELEMENT PROGRAM FENA. 3D 


5.1 Test Problems 

In order to make some assessment of Krishnayya's 
program FENA.3D, two types of test problem were investigated: 

1. A perfectly regular homogeneous dam, to help 
evaluate the capabilities of FENA.3D. 

2. A thick cadlumn comprising some “odd" elements, 


to evaluate some aspects of element performance. 


A wedge-shaped dam in a V-shaped valley: Data and 
mesh as described by Lefebvre Duncan, and Wilson (1973) 
were prepared and run in order to gain some experience with 
the program FENA.3D and to compare some results with a 
published work: The structure is shown “In Fig. 5.1. A 
coarser analysis (4 lifts) was also performed. Useful 
information to assist in prediction of the computation time 
and total machine costs for further analyses was also 
obtained from this study. 

Results for the 8-lift analysis are compared with 
the results of Lefebvre et al. in Appendix A. Some discrepan- 
cies occur, and may have arisen from several sources. For 
a further discussion of the matter the reader is directed to 
the notes in Appendix A. The writer concluded that with normal 
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care in the preparation of data the results from FENA.3D 


would be quite satisfactory. 


A simple thick column? A series of three small 


analyses were run on the simple steel column shown in Fig. 
5.2. In theory the vertical stress should be constant at 
4.444 Ksi but the effects of non-conforming elements and 
node-ordering of the elements was investigated. Such effects 
arose during construction of the Mica mesh when it was found 
impossible to avoid non-conforming surfaces (for example a 
hyperbolic paraboloid versus two triangular faces in Fig. 5.3). 
This physical problem has a numerical counterpart in element 
stiffness formulation, and even though all such situations 

in the Mica mesh occur in the least interesting areas, it 

was felt necessary to place some bounds on the inaccuracies 
which might accrue. 

Additionally, the effect of the nodes of an element on 
the numerical integration procedures was examined for condensed 
elements. Condensed elements have one or more nodes coalesced 
as shown in Fig. 5.4. Though a basic order must be obeyed, 
the starting position can vary giving different repetition 
possibilities in numbering of these elements. 

Results of this study are shown in Fig. 5.5 and errors 
of up to 22% occur. These errors are not systematic in any 
simple fashion and, while serious enough, the positions in 


the dam where the situation occurs are: 
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FIGURE 5.2 SIMPLE THICK STEEL COLUMN 
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(a) not important from the viewpoint of the interest 
in stresses, 
(b) liable to certain inaccuracies in any case due 
to artificially stiff boundary zones, 
(c) unlikely to have significant influence on the 


structure as a whole, because of limited nodal connectivity. 


Conclusions! from test problems: The tests clearly 


indicate that in the Mica mesh, where elements are always 
skewed and often highly distorted and some are non-conforming 
and/or condensed, some numerical inaccuracies must be 
expected and all results accepted with certain reservations 
as to accuracy. The effects of numerical roundoff in the 
computer must be considered realistically too, the IBM 360/67 
has only a 32-bit word length and all computation was 

done in single precision. However, consistent and reliable 
results will be obtained in the areas of chief interest, 
namely the central transverse section and the deep zones of 


the core. 


5.2 Homogeneous Embankments 


General: Kulhawy et al (1969) and Clough and Woodward 
(1967) have made very significant studies of homogeneous embank- 
ments as an aid to furthering the engineer's understanding of 
the methods and processes, and alsotogive a feeling for the 


numbers involved. 
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It was felt appropriate to run a series of 

analyses assuming that Mica Dam was homogeneous, to give 
a simple property (easily stored in one's head) which 
roughlysexhibits thesfield deformations: « Thus; a "shell" 
property was determined which matched settlements in the 
shells of the maximum transverse section (at Station 22+50), 
and a "core" property by matching settlements in the 


longitudinal core section. 


pone Lit property of Mica Dam: The linear elastic 


parameters chosen for zone M2 were used ina 1 lift 
("gravity turn on") and 5 lift analysis. The displacement 
pattern of the 5 lift matched field settlements well, and 
correlation factors for the 1 lift analysis were determined. 
Using similitude (with approximations) all subsequent 
homogeneous analyses were made in 1 lift, and the settlements 
converted to the 5 lift pattern. Variation of Poisson's 
ratio had a small effect on settlements and maximum principal 
stresses and a larger effect on horizontal movements and 
minor principal stresses, as expected from the works cited 
above. 

To match settlements in the shell as well as possible, 


the following parameters were used: 


D 9800 Ksf U = Diez a 


3000 Ksf 


i) 
II 


or K 5810 Ksf 
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‘Core" property of Mica Dam: It was found that core 


material properties derived from laboratory oedometer tests 
or field MV gauge data did not match but considerably over- 
estimated core settlements. A homogeneous embankment which 
matched settlements in the core as well as possible had the 


following parameters: 


D 


4900 Ksf u = O35 


oha K 3400 Ksf G 1134 Ksf 


The discrepancies between these values and derived values 
(Section 4.5) arose because of the homogeneous structure. 
Further discussion of these results may be found in Section 


Gail 


5.3 Analysis of the Zoned Embankment 


silupleylinear analysis: A 5-lift analysis; assuming 


rigid boundaries, wasS run using seven materials with 
properties as described in Section 4.5. Very good general 
agreement was obtained with the field deformation pattern. 
Relatively, core settlements were underestimated and shell 
settlements overestimated. Stresses in the core were so 
low that it was decided to isolate the effect of arching 
due to steep valley geometry if possible by using a plane 


strain analysis. 
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Plane strain analysis: The mesh between Stations 
22+50 and 25+50 was isolated in plane strain by fixing all 
nodes on or between these sections against cross-valley 
movement. Surprisingly enough, the same linear elastic 
parameters as used above gave results very similar to the 
simple 3D linear analysis. Settlements were generally 


slightly overestimated. 


Linear analysis with bedrock movements: Using the 


same mesh and properties as for the simple linear analysis 
bedrock movement during construction was simulated by applying 
field-measured increments at boundary nodes on bedrock, as 
discussed in Section 4.6. There was little change in stress 
distribution but a slight improvement in matching field 


settlements in the fill. 


Multilinear analysis including bedrock movements: 


Using the "average stress" approach, allowing the deformation 
moduli in the elements to vary with stress state (as 
described in Section 4.5), and applying the bedrock 
movements, another analysis was made to determine in a simple 
fashion the effects of increased oobi ican iaa of material 
behaviour. In general stress distribution was not drastically 
affected and settlements mostly agreeable with field measure- 
ments. 

All the results processed from these analyses, and 
some conclusions which can be drawn, are discussed in 


successive chapters of this thesis. 
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CHAPTER VI 


RESULTS OF THE 3D MICA DAM ANALYSES 


6.1 The Homogeneous Embankment Analyses 


As discussed in Section 5.2; the results of the 
homogeneous embankment provide a guide to the sort of 
phenomena which occurred as Mica Dam was constructed. Zones 
of tensile minor principal stress were observed in the areas 
expected (Krishnayya, 1973a); Covarrubias, 1969). These 
were of limited extent and the coarseness of the mesh inhibits 
any more than passing interest in this. A 5 lift analysis of 
the homogeneous embankment did not display any tensile 
stresses. Examination of the displacements and normal 
strains in the longitudinal direction at transverse Station 
22+50 definitely ruled out plane strain behaviour at this 
area, which would be the assumed plane strain section for 
2D analyses. Regardless of this, such movements were generally 
about an order of magnitude less than the corresponding 
Maximum movements and strains at this section. It is there- 
fore concluded that, relatively speaking, an assumption of 
plane strain behaviour at this section is valid for 2D work, 
supporting the argument of Skermer (1973). Stresses calculated 
for the homogeneous dam are discussed underneath, and the 
reader! is) referred back to Section 5.2 for a discussion of 


parameters used. 
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6.2 The Zoned Embankment Analyses 

Figures 6.1 to 6.4 show the agreement of computed 
settlements and field measurements for pertinent MV gauges. 

A field settlement plot consists of lower bound and upper 
bound lines. The lower bound is the cumulative settlement 
during construction periods only, while the upper bound 

is the cumulative settlement for construction and winter 
periods. The smaller indicates "elastic" response and the 
larger includes creep and consolidation effects. Plotted on 
the field curves (which are simplified for case of inter- 
pretation) are settlements for the cases of linear elastic 
3D and multilinear elastic 3D (both with bedrock movements), 
simple elastic 3D and simple elastic plane strain 2D. 

Since some attempt was made to include consolidation 
effects in the core in the modulus selection process, one 
would hope to see calculated core settlements tending towards 
the upper bound of field measurements. Since elastic 
parameters do not attempt to include creep effects, calculated 
shell settlements would preferably identify with lower bound 
field measurements. In general, the overall pattern of 
calculated settlement agrees well with field measurement. It 
could be argued that lack of detail in the mesh restricts core 
settlements and helps explain the general underestimation of 
field values in this zone. However, the multilinear analysis 
considerably overestimates core settlements in the softer 


(1969) material, and the explanation advanced is that 
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FIGURE 6.1 SETTLEMENT PROFILES: MV7 AND MV8, STA.22+50 
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FIGURE 6.2 SETTLEMENT PROFILE: MV9, STA.22+50 
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setHement profiles, simplified from field data. 


FIGURE 6.3 SETTLEMENT PROFILES: MV'S 10, 11, 12; STA.22+50 
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FIGURE 6.4 SETTLEMENT PROFILES: MV5 AND MV15, STA.25+50 
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87 
insufficient experience and perhaps not the best selection of 
test data caused this. Also, since shear failure was 
considered on an element basis (as described in Section 3.3, 
using the Mohr-Coulomb criterion) and high shear strengths 
are mobilized chiefly in the core-shell interface zone the 
finite element representation may well be too approximate. 
Perhaps shear failure should not have been considered in order 
to gain a more representative analysis with such a coarse 
mesh. Shell settlements were generally overestimated, and 
this may perhaps be blamed upon the field deformation data, 
which could not be separated into constructional and total 
strain increments. In general the settlement results are 
good and since settlements are good we should expect that 
stresses, which are less sensitive to changes in moduli alone 
(Clough and Woodward, 1967), should be quite realistic as 
compared to the structure. 

The plane strain 2D settlements overestimate the 
3D settlements from an analysis using identical element 
subdivision and material properties. This is important as 
it shows the effect of 3D geometry which allows deformation 
modes not allowed in 2D plane strain. The 2D plane strain 
analysis used a section of the 3D mesh. Palmerton and 
Lefebvre (1972) indicate that comparable 2D plane strain 
analyses using 2D and 3D elements show less settlement in 
the 2D element case. Thus a true 2D plane strain analysis 


which gave close agreement with field settlements may be 


eo) SD ee, See eee 


en teed sii gon aqeAasy 
aew Stukhial weere esate | oBtA | 
6.€ nobjee2 st Betieeesh en) abasd snemahal as mo 
aifunerte wade pid bas (nodwaditD dot fuod-1tom ott eri 
add eno= sostrsint Llede-o'%00 afd ak yiterds bes idom S38 
.etboltxoiqus ont sd Iiow yem jos tecnsae7ge2 snemste eatnlt 
webro Hi Ss4hbRanoo. aesd even 400 pivots Suufist xsefe eqediet © 
327509 5 Hove diw Heetbns oui seinentgen szom & nisp oF 
bis ,Secamitgausvo ylisiensp eIpw edhemel aoe {lesa . stam 
,646b aclipmoleh Bieiz 4c) cogt, emake ? egedisq yem eid? ] 
lates bns fenotinuypaenoo ost: fia teysqee sd ton biwoo dod 
ei6 wilgess inemaligve df3 fexsnoe of .eonemexont aiage 7 
tern tosqxs oivoda sw Soop ets bsromeiaem sonia baie Boo 
enole ifvubom ni aspastio 92 evidiense esaf ozs dotdiw reoneoise 
es sitctise1, stiup od Bluofa . (Teel enesciaisamtel brs 
vetnemmaal os 
ott sismitestovo evnemeitsee OS abate -enplg edt 
tnemeais [estinebi priav eieylens 16 aor? es 
as Jn&3rtoqmi'ei 2int seers 7zeqesq ‘sizesam bas no. 


| atsvse enefaq. OS edt .nisade onslq GS at bewolis Joa 
bee aosasmiad .fesm 06 ef? 20 sokdoew & Beau 
nisizs saeig OS edetsqmoo tedt eraolfint (svtet) 

ni Jnempisde® s2o! wore etnemole Gf Bas GS eatev assy 

siaylans o, ensiy GS oust 6 gudT 9889 anomeis aS. 


7 
ps Out 


HORIZON TAL, MOVEMENTS 


= 


\e 


line of 


= 


7 


MESH SCALE 
FEET 


3 

ae 2 

Els Ee \s 

Sy = ul bf 

cle , Ww e: 

ie) j Si 

+} io 

Gg 

33 ae 

Us > < oy 
oO) > 

rad a 2” 

Ole 

| 5 : 

al 

fe es] \ . 

_:| 3 XN 

uJ ie 2 

ale =~ 


= 


FIGURE 6.5 PATTERN OF HORIZONTAL MOVEMENTS WITHIN DAM 


Lui’ DY 


BROMGSOMLYT WOAEAEALS AI 


bYLLESBY GE 


RTCGWE &*. 


89 
quite representative of the field behaviour. Mica Dam 
was expected to demonstrate obvious 3D effects due to the 
changes in valley direction and slopes. That these results 
are not very evident may be a measure of the coarseness of 
the 3D analysis. Also, it may demonstrate the principle of 
St. Venant, the dam is a huge structure and influences from 
a local area are not felt at larger distances. 

Horizontal movements measured in the dam are very 
small, and no direct comparison can be made with computed 
values. This is unfortunate since the role of Poisson's 
ratio is quite important. A pattern of horizontal movements 
atrapprox. el .@2150%.rt. (the\end of construction in 1970) 
is shown in Fig. 6.5 to demonstrate typical horizontal 
movements expected after the embankment was complete. The 
lack of confirmation of the chosen values of Poisson's 
ratio restricts the confidence which can be placed on 
computed values of minor principal stresses, particularly 
in the core. However, the expected variation in stress 
would not likely be more than 10-20% of calculated values, 
in the writers' opinion, based upon assessments of Poisson's 


ratio variation in the homogeneous embankment analyses. 


6.3° Stress Distributions 


Core: Figure 6.6 shows a plot of element stresses for 
a column of core material between stations 22+50 and 25+50. 


There are only two elements spanning the width of the core 


20 aighoatiy oft ecamanchoh 
mor? a ors exutanyt®: 
gsonndeih seprel Je 

yrev S36 mab sda ni. hace bale 


besuqmos gditw obam 4d. 463. 


a'nnseted 3d olex sit conte sdeauaze%no ni iat 
asremevon iasgasbrad Qo su7eg 4 -isiogad ani 
(OFHL ns aoitou7 tenes an boe:adt) 23 ovis ke xotaanite 
fasnosixed Lewray? Sap payee 39) 2.8 sna ha 
ait .asialgmoo 25w Tagminsclns ons Sealed 
s'nowatot- So seulev st ony 6 

no bevsle sd mso ASEAW « Zi 
viselporsaag iaabbenianona 
neortte Al nol deinen ae anes 
eoulsv bolte tastes ta ‘Hoses, nails panei 


2'aceatot Io etnamedanes nero aed - aa 1. 
-aoay inns sasminsdmes evoonseumet! ef zo 7 
| | | 3; is 


- 
a 


1. 


, kd 


MAX. 


PRINCIPAL STRESS (Ksf) st 
2500 _1° 20 40 60 80 100 
ELEV. \ \\ 
(ft) A Upstream and Downstream Elements, 
N Core Column: Average of 
2400 
\ 
\ 


Centroid Stresses 


bly 
2300 


lineav, incl. bedrock mvts. 
multilinear, incl. bedr. mvis. 
linear, 2D plane strain 


homogenews dam, 
2200 


Skermer's 2D analysis 
(1974) 
2100 


2000 
1900 


1800 


FIGURE 6.6 


STRESSES IN THE CENTRAL CORE SECTION 


OL 
from upstream to downstream shells, and the displacement 
pattern permitted is not fully representative of field 
behaviour. Thus it is stated clearly at this stage that 
the writer considers all stress results to represent average 
effects only, or a good first-order approximation at best. 

The element major principal stresses and vertical 
stresses are close ‘in value and for clarity only the 
principal stresses are plotted. The indication is that 
the prime mode of deformation is vertical settlement 
restrained by the mobilization of high shear strength along 
the core-shell interfaces, particularly the downstream side. 
This verifies an intuitive assessment. For the homogeneous 
embankment under various combinations of moduli, essentially 
the same vertical stress variation as that shown in the 
figure resulted. ‘The conclusion is that no appreciable 
stress transfer occurs to about elevation 2150 ft, below 
which valley arching relieves about 30% of overburden stress 
at elevation 1860 ft. (approximately the old river level). 
Unfortunately a plane strain 2D analysis was not carried 
out using the homogeneous parameters, although the writer 
believes that stresses would have been very similar to those 
mentioned above. 

In the zoned fill, however, the prime factor 
determining stresses in the core is the change in material 
stiffness from core to shell. Even in the upper reaches of 


the core stress transfer is about 40% of overburden; this 
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92 
value really starts to increase below elevation 2150 ft. 
with the influence of lateral (valley?) restraint, and 
reaches 65% at elevation 1860 ft. Somewhat surprisingly, 
stresses from the plane strain 2D analysis agrees very closely 
with 3D results. This indicates that valley restraint is 
minimal and that lateral restraint has only a secondary 
effect on stress distributions in this case. 

The really significant aspect is that Mica Dam is 
admittedly a 3D structure, and theoretical work of previous 
investigators (Palmerton and Lefebvre, 1972; Lefebvre, Duncan, 
and Wilson, 1973) indicates that 3D effects may be quite 
important as regards stress distribution. The analyses 
described herein imply that although 3D effects may be present, 
the governing deformational behaviour can be represented 
satisfactorily using a 2D plane strain analysis. Thus, 
increased detail and the sophistication of stress-strain 
response which can be incorporated into 2D analyses should 
yield satisfactory results for engineering purposes, and only 
for a study of cracking potential would the 3D analysis be 
required. 

A completely different, alternative conclusion is that 
the 3D analyses described are too coarse to provide any 
representative results. This conclusion is rejected by the 
writer until satisfactorily proven. For comparision, Skermer's 
(1974) stress results using a 2D analysis are also shown on 
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Shell: The stresses computed for the shell are 
not much different from what one might expect. The upstream 
shell essentially sits tightly in place over the core, and 
the core widening at the base causes high stress levels as 
the shell overburden is thrown into a smaller bearing area. 
The downstream shell essentially sits tightly in place, 
partially under the core, and participates more strongly in 
core-shell interaction than the upstream shell. The weaker 
zone M2DI exhibits no appreciable stress transfer at depth 
because it holds the burden of core-shell interaction, 
and also because of the wedging action of the foundation 
excavation. The shell stresses for the linear analysis 
including bedrock movements are shown in Figs. 6.7 to 6.10. 
It should be noted that these stresses are quite comparable 
with Skermer's results, although he did not model material 


; 1 
differences to the degree adopted in the 3D work. 


Overall stress distribution: For the sake of 
comparative evaluation stress distributions in transverse 
and longitudinal sections are given in some detail for 
various analyses. They are all grouped together for convenience 
in Appendix D. Data is plotted using engineering judgement 
and without making any allowance for artificalhly stiff boundary 
effects. Because of the coarseness of the mesh there is 
considerable scope for personal interpretation of contour 


patterns, although the writer has endeavoured to be consistent. 
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Element stresses are calculated at the element centroid and 
nodal averaging procedures were not used because averaging 
the values of such a coarse mesh was found to be unsatis- 
factory. Stresses in the transverse section shown as 
Station 22+50 are actually in the elements between Stations 
22+50 and 25+50. Stresses in the longitudinal core section 
are averaged from the upstream and downstream core elements. 

Attention is drawn to the low stresses in the core, 
a serious situation which may lead to hydraulic fracture. 
This is described further in Section 7.2. As an estimate 
the writer believes that stresses at some point (particularly 
in the core) would not vary by more than 10 to 20% from the 
results of his anslyses. The stress situation in the core 
of~the dam'is not as.satisfactory as it could be, that is 
the simple truth of the matter. Whether there is any 
practical significance in this point is a matter for ongoing 


research and debate. 


6.4 Stress Paths of Elements 

The elastic parameters that were used to represent 
deformational behaviour at varying stress levels only 
crudely represent the physical processes of the material. 
Some insight into the approximations dealt with can be found 
by plotting an element stress-strain path on a background 
of (say) oedometer tests for that material. Also, the plot 


of octahedral normal stress versus maximum principal stress 
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May yield some insight into the behaviour represented by 
Poisson's ratio. Plots for a typical core element are shown 
against oedometer data for the core material, for the linear 
analysis which gave the best agreement with field deformations, 
on Fig. 6.11. On Fig. 6.12 are octahedral stress/principal 
stress plots for several core elements. Rough correspondence 
of moduli are shown, and theeffective value of Poisson's ratio 
is close to constant at the assumed value. The same is 
basically true of shell zone M2 although the effective value 
of Poisson's ratio drops off a little with increased stress 
level. (Figs. 6.13 fand 6.14) . 

The conclusions which can be drawn from this are 
limited. One would expect that the material might exhibit 
a wider variation in Poisson's ratio, but that some average 
value may be broadly representative. One would also expect 
that materials might exhibit considerable variation in 
moduli, but that some average modulus could describe this 
adequately. It is not absolute values of elastic parameters 
which govern stress distribution but rather the ratio of their 


values. 


6.5 Comments on Some Important Influential Factors 


Shear failure of elements: In accordance with the 
work of Kulhawy et al (1969), elements in which computed 
shear stresses exceed the calculated failure shear strength 


(using the Mohr - Coulomb failure criterion) are assigned 
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Test Data from T.C. Law (1974) | 
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we Bost Fit to MV Gauge Data : Shell Zone M2 


(1969 ~1972) 
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FIGURE 6.13 STRESS PATH: SHELL ZONE M2 ELEMENTS 
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low values of shear modulus for subsequent stiffness 
calculations. In principle this feature is worthy of 
inclusion, but there is some doubt in the writer's mind as 
to its suitability in a coarse mesh, as Krishnayya (1973a) 
showed that failed zones are governed to some extent by mesh 
size. In these analyses too much material may have been 
assigned "failed" properties, although preliminary estimates 
indicated that the "failure - no failure" option would have 
less than 10% influence on stress distribution. Elements 
which failed in shear were either highly distorted elements 
in poorer and unimportant areas of the mesh, or elements 
at abrupt changes of slope in the longitudinal core section. 

Ideally an analysis which considers shear failure and 
which also uses the "no tension" approach (Zienkiewicz et al, 
1968) would be desirable. The writer doubts whether either 


is practical or realistic in the present analyses. 


Sophistication of stress-strain relations required: 


The evident success of the finite element method as applied 

to fill construction analysis strongly suggests that, for 

the stress paths considered, representative behaviour can be 
obtained by simple approaches to nonlinearity, or representative 
linear parameters which have to be chosen very carefully. This 
is an encouraging result for it makes the method much easier 

to use in practice, and allows efforts to be concentrated 


upon proper simplifications for future research. 
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No one should ever believe the results of a finite 
element analysis verbatim without careful assessment of the 
considerations which went into choice of parameters. This 
work on Mica Dam reinforces the opinion that finite element 
analysis is now a useful design tool in geotechnical 


engineering, when used wisely. 


General accuracy required and claimed: Inaccuracies 


of many types were bound to accrue in the Mica Analyses. 
Coarseness of the mesh, element performance, representative 
stress-strain behaviour, options for shear failure, and 
numerical roundoff all are present. No serious probabilistic 
approach to assessment of overall errors is possible. The 
maximum likely error is estimated as about 15-20% and it is 
claimed that it would not be justifiable to seek a reduction 
of this figure. Furthermore, these results are quite 
acceptable for practical soils engineering purposes. To seek 
greater accuracy would be unrealistic and contrary to well 
established practices. Nothing should ever supplant the 
observational method as the mainstay of soil engineering 


practice? 


Cost of 3D analyses. Some remarks concerning the 
costs involved in making the analyses may be found in Appendix 
E. 

1 oxermer (1974) used only three linear elastic 


materials in a rather arbitrary fashion, and his results must 
be treated with caution. 
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CHAPTERS Vit 
SOME IMPLICATIONS OF THE STRESS ANALYSES 


The end of construction showed that the Dam was 
performing satisfactorily. Embankment movements were 
smaller than anticipated and the structure was deemed 
successful. One of the consequences of these small movements, 
perhaps, is the exceptionally low stresses in the core due 
to the high degree of stress transfer between core and 
shells. In the light of these low stresses the following 
questions arise: 

(1) What effect will the reservoir have, on first 
filling and in the long term? 

(2) What possibility is there of hydraulic fracturing 
of the core? 

(3) Do these low stresses reflect any concern in 
the safety of the structure with respect to earthquake 


loadings? 


foie tibat er wllingeols the Reservoir 


Nobari and Duncan (1972) undertook a detailed 
investigation of the effects of reservoir filling upon an 
embankment. Complex movements had been observed in many 
structures, causing cracks due to differential settlement and 
hydraulic fracturing. The basic mechanisms are: 
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(a) compression due to wetting of the foundation, 
causing settlement which is non-uniform due to the progress 
of the wetting front, 

(b) compression of the wetted portion of a 
homogeneous or zoned fill, causing settlements and upstream 
movements and spreading possibly, 

(c) downstream movements, due to the water load 
which is greatest during the later stages of filling. 
Embankments usually exhibit individualistic variations of 
the above based upon local conditions and materials. 

Two basic mechanisms take place as the reservoir assumes 
design level. Firstly, before significant seepage patterns 
are developed, the water load will increase the total stress 
state in the core. Depending upon the existing construction 
pore pressures and the core material response (say B) the 
effective stress state in the core will usually increase 
somewhat as well. Secondly, as the steady state seepage 
pattern forms, there will be some changes in effective 
stresses in the part affected by seepage. These two 
mechanisms are not likely to be easily separated or identified 
in general because of any redistribution of stresses during 
wetting and because of the masking effect of construction 
pore pressure dissipation. 

Mica reservoir (Lake Kinbasket) started filling in 
April 1973 and reached a level well below minimum pool before 


discharge requirements under the Columbia Treaty reduced it 


‘iden 


ase ttage bts atremtjses pateuso villi benot to 
bEQL 1970W oft of esh .ednemevom sseczteawob. (>) - +) =, 
Me LE LS Do coypde s926d ade puixeb tesdsexp ab diohiw 
te atotasivtey site isibevionk tididve yiisuen esasmiasded 
-2ipiig¢am bas aneitibaes fmeol nega beasd evods oft 
aemuzes tLoytade 204 ge acnlq si6o eteinsdoom oliasd owt bd hoa 
enssitsq spedsee tnsobilapie sueted ,ylse1it -level-npiesb 
seais2 Leto? ofY sanexoni Lily Baol astew eddy, Deqoiere® sis 
soijoutianos poitelxe es nego gnabasged «s100 Sav at o3sde 
oft @ yeu? sanomser Isitedem st0° sft has aetueestq 230q 
saneazoni ylisveu ILiw e160 edt ni stase seo1se svisseite 
opaqase stsde ybedse ofd we .yibnovea + .Llew as Jesivomon 
eVisostis mi cant eine ad Ly oxadd senso} meding 
ows sasiT .epmqses yd besootis 3steq eft nt soesoxda 
hoititashi 16 bedexeqss yliese ed ot yledil — os 
prbtad esagotie to avi tual atathes vie Yo ‘seusoed 


ih 


noisovssenae 20 sonia entaeen ods te saan ee 


108 


STRESS 


excess porewater pressures 


construction reservoir fi lling 
completed. commenced. 


FIGURE 7,1 DIAGRAMMATIC RESPONSE OF EARTH PRESSURE CELL 
TO RESERVOIR FILLING 
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fartherss It may filletorwithins50 feet of maximum pool by 
mid-summer of 1974. Preliminary results from some earth- 
pressure cells monitoring horizontal stresses in the core 
responded as discussed above and shown diagramatically in 
Bag ime: (Skermer, personal communication, 1973). It is 
Loosearmly;eatethedstage of writing; tondiscusstreservoir 
PvEPing at Mica Dam in further detail. The subject requires 
a separate research program due to the complex material 
behaviour required for finite element modelling. A good 
outline of the approach required is that given by Nobari and 
Duncan (1972). 

The special care taken in design of the Dam should 
ensure that reservoir loading will have a beneficial effect 
on the stress state in the core. Many embankments have been 
constructed embodying the arch principle in the core but 
to the writer's knowledge no specific study or quantitative 
assessment has ever been made of curvature, If for various 
reasons some cracks do form allowing the percolation of water 
the core material is "self healing”. The initial erosion 
process washes in materials which will be caught in and 


eventually seal up the cracks. 


7.2 the Likelihood of Hydraulic Fracture 


Blight (1973) presented a means of assessing the 
stress transfer in a narrow core or trench. Some of the 


parameters required are obtained by "reasonable engineering 
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estimates". One important parameter is the coefficient of 
lateral earth pressure. Blight makes use of some high values 
Oia 1 5e0Gn2. Wilk Kis calculated using elastic theory and 


Poisson's ratio, much lower values may be obtained. Making 
estimates of the parameters the following results were 


obtained for the centre of the Mica Dam core at elevation 


36 Oeft2 
Ky equiv. i Stress transfer (3) 
eae) 0.66 76.0 
10 Oe 57.0 
fa 0.33 3 70 


Blight's approach was not successfully applied because 
it cannot agree with the finite element calculations. The 
estimates may not have been satisfactory, of course. 

Blight cited the example of Balderhead Dam (Vaughan 
et al, 1970) where a hydraulic fracturing phenomenon was 
observed. His results indicated that piping failure at 
Balderhead was initiated by a fracture in the core which 
occurred during construction. This conclusion is not, shared 
by the original reporters nor by Sherard (1973) in his 
summary of embankment dam cracking, so that Blight's 
analysis must be viewed circumspectly. 

Skermer (1974) reported the results of water inflow 
tests in the MV 15 gauge at Mica Dam. His evidence supports 


the reliability of his finite element analysis which yielded 
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very low stresses in the core. Data and water pressures 
are taken from his work, and minor principal stresses from 
thes;present 3sDanalyses, and shown on Fig. 7.2... Reservoir 
level was at el. 2218 ft. Water Wevel inthe MVegauge could 
pot bevraised above et. 2185 fe. despite copious mnrlow. 
Undoubtedly hydraulic fracturing occurred at some point around 
the casing, most likely below el. 1950 ft. where computed 
minor principal total stresses are less than water pressure. 
Alteratively the water could have escaped through bedrock 
as fissures opened up under the high water pressure in the 
gauge. To the writers' knowledge no information exists as 
to whether the reservoir level itself caused any hydraulic 
fracturing. According to the above studies it could have. 
Perhaps the self-healing properties of the core have closed 
any such fractures. The questions regarding hydraulic 
fracturing of the core of Mica Dam are not fully answered at 
this stage. Under the circumstances, though, leakage due to 
hydraulic fracturing is not considered likely. 

The phenomenon of hydraulic fracturing is excellently 
covered by Sherard (1973). Deliberate fracturing in boreholes 
has been utilized in the high pressure Soletanche "sleeve 
grouting" procedure for grouting alluvium under dams, while 
Bjerrum and Andersen (1972), Vaughan (1972), and Penman (1972) 
have discussed the use of induced hydraulic fracturing to 
estimate the earth pressures in materials (in a similar 


manner to that illustrated at Mica Dam). Sherard observed 
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FIGURE 7.2 INFLOW TEST IN GAUGE MV15 
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that the use of fracturing in boreholes, particularly involving 
large water inflows, may lead to a dangerous situation where 
none previously existed. He pointed out the different 
stress situations of borehole-induced and reservoir-induced 
hydraulicetiracturing. 

Sherard, et al (1963) present comprehensive data 
gathered for earth dams in North America, which shows that, 
after the non-geotechnical problem of flood-overtopping, 
piping failures comprise the greatest proportion of unsatis- 
factory earth dam performance. Evidence from the World Register 
of Dams supports this conclusion. Sherard, Decker, and 
Ryker (1972) described the piping failures of several low 
embankments and concluded that hydraulic fracturing initiated 
the fissuring which led to the piping. Evidence from finite 
element analyses supports their argument. It may well be that 
many more of these "piping failures" are introduced by 
hydraulic fracturing mechanisms, unrecognized because of the 
limited knowledge available concerning stress states in 
embankments. Very few confirmed cases, where hydraulic 
fracturing has initiated leakage in large dams, have been 
reported and discussed in a thorough fashion. Finite element 
analyses may play a fundamental role in furthering the 


understanding of this process. 


7.3 Earthquake Design of Mica Dam 
As discussed by Sherard (1967) and Seed (1973), there 


is a serious lack of information concerning the earthquake- 
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resisting design of high embankments. The limitations of 
analytical procedures are well illustrated by Seed; who in 
summarizing states that "Previous performance of earth dams 
during earthquakes gives no grounds for complacency regarding 
the ability of well-constructed embankments to withstand 
the effects of a major earthquake, although, if suitable 
design measures are taken, they give no ground for concern." 
Both authors state that the most effective features which 
can be employed are: 

(a) a core of sufficient plasticity and of good 
self-healing properties, 

(b) good filtering properties at the core-shell 
interfaces, 

(c) provision of adequate freeboard, crest width and 
erosion resistance against overtopping, 

(d) minimising of the risk of liquifaction of the 
upstream shell. 
Sherard also gives a summary of the best properties for core 
Material based upon experience and assessment. 

Intuitively the most satisfactory stress distribution 
in an embankment is as uniform across discontinuous zones 
as possible, and this can be achieved with core and shells 
of equivalent stiffness. This is invariably impossible to 
achieve in practice because the upstream shell is required 
to be as stiff as possible to minimize the risk of liquifaction 


(Casagrande, 1965). At Mica Dam this was done, and reliance 
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placed upon the self-healing properties of the till core and 
the good filter-transition properties of the sandy gravel 
downstream shell to safeguard the core against any crack 
formation (Skermer, 1974). 

The effects which an earthquake would have on Mica 
Dam, bearing in mind the stress state in the core, are not 
Clearly understood by the writer. With the provisions that 
have been made in design it may be assumed that the structure 
would be as safe as it can be regarding reasonable design and 
extensive experience as against foreseen circumstances. 
However, the 3D finite element analyses which have been 
reported in this thesis indicate thatwuah2D analysis should 
be quite representative. It would be possible to undertake 
a 2D dynamic response analysis, given information concerning 
a possible earthquake, to determine after some fashion stress 
or strain data which may prove valuable in investigating the 
possible effects of a shock on the structure. The knowledge 
of stress and strain data under dynamic loading would be 
even less substantial than the static case (which is quite 
limited) so that properties used in such analyses would need 
to be chosen with care and results treated with caution. 


Such an exercise may be quite instructive. 
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CHAPTER VIII 
CONCLUSIONS 


Four general classes of conclusions can be drawn 
from the 3D finite element analyses described in this thesis: 

(1) Scope, implications, and acceptability of this 
work within the range of finite element analyses which can 
be used at Mica Dam. 

(2) The role of these analyses in discussing the 
field behaviour of Mica Dam. 

(3) The general validity of such an approach and the 
sophistication required for analysis. 


(4) Specific recommendations for future work. 


Scope, Implications and Acceptability 


It has been shown that the coarse 3D analysis, using 
stress-strain data derived mostly from the field, has 
successfully reproduced aspects of field behaviour. This 
success of the model suggests that the stresses calculated are 
representative of conditions in the embankment, and some 
field evidence may confirm the stress results. As discussed 
in a previous investigation (Eisenstein, Krishnayya, and 
Morgenstern, 1972b), there exists useful information which 
can be utilized for assessing the cracking potential of the 


core of Mica Dam. 
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The analyses suggest that despite the markedly 3D 
aspect of the structure, the constructional deformation and 
stress behaviour can more economically and just as success- 
fully be obtained from 2D plane strain analyses. 

For the materials and loadings considered further 
evidence has been presented that gross but thoughtful 
simpliciations in material behaviour can be made and still 
produce adequate information for design and monitoring purposes. 

Reasonable inferences as to the effects of reservoir 


filling can be made from the present results. 


Role) of) theyAnalyses,in FieldsStudies 


Since most deformational data was obtained from 
field measurements the agreement of field and analytical 
settlements is a mark of self-consistency of the data only 
and indicates that in the core, where laboratory data was 
used, the finite element results vindicate the simple deform- 
ational moduli selected. 

If leakage due to hydraulic fracturing is detected, 
as seems highly unlikely, the analyses should provide very use- 
ful information. In advance of any further studies, the 
results will shed some light on any unexpected behaviour 
observed during reservoir filling and possibly earthquake 
loading. 

Some information is available concerning stresses 


and movements in areas where conventional 2D analyses may not 
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be justified, such as zones of potential cracking. It is 
important to recognize the coarseness of the present 3D 


coverage in this respect. 


General) Validity of the Approach 


Results of these finite element analyses indicate 
that even with such gross simplifications in stress-strain 
representation and mesh geometry successful finite element 
modelling of Mica Dam occurred. The lack of detail should 
always be borne in mind, but also bearing in mind the 
accumulation of finite-element-analysed case histories now 
available, good average data which can be applied in 
engineering design and evaluation are now available. 

Apparently, the simple linear elastic incremental 
analyses can provide satisfactory modelling of the structure 
and for general design purposes any further sophistication 
is unwarranted except in some severe cases. The data which 
would be available before construction may cause serious 
disagreements between predicted and measured behaviour, 
throwing light on the validity of sampling and testing 
procedures rather than the calculations themselves. 
Obviously there is much scope for effective research in 
developing reliable test methods to represent constructed soil 
behaviour. 

It is worthwhile to bear in mind the comment of 
Smith and Kay (1971) concerning the apparent success of very 


crude approximations to true soil behaviour. When applying 
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the finite element method to phenomena where the basic 
deformational response is pseudo-elastic, we can expect 
reasonable results using a pseudo-elastic analysis. The 
effects of consolidation, creep, yield, hysterisis and 
dynamic loading are a whole world of complexity different, 
and at the present stage cannot be realistically employed 
with great reliability in finite element analyses for design 
purposes. Parametric studies in these fields will continue 


to be of great use, however, as research develops. 


Specific Recommendations for Future Work 


The effects of reservoir filling will likely be of a 
3D nature Owing to valley geometry and core curvature. Such 
an analysis, taking into account the wetting of the upstream 
shell, will be valuable and will be 4 timely contribution 
to the state-of-the-art not just regarding the mechanics 
of wetting and filling but particularly in assessment of the 
effects of curvature and other procedures to optimise wedging 
action and minimise serious crack formation. Also, due to 
the great height of the embankment, detailed finite element 
analyses in 2D using the best means to account for 
constructional behaviour with test data will contribute 
significantly to the general confidence in and development 
of representative test procedures. 

Earlier work (e.g. Palmerton and Lefebvre, 1972) has 
indicated that in certain situations a 3D analysis would be 


economical enough and useful enough in assessment of 2D work, 
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because of the extra information and more representative 
modelling. The restrictions on 3D analyses are the data 
preparation and computer resources which must be involved, 
and no quick easing of these restrictions can be envisaged. 

Finally, it is readily apparent that the finite 
element is a useful and practical design tool when used 
sensibly within the limits of its capabilities. Future 
researchers should bear in mind always the questions of 
practical application, otherwise vast amounts of computer 
resources and energy supply are being expended in playing 
with sophisticated toys. Users of the method in industry 
must always educate themselves to the developments which 
occur and be prepared to change their concepts of how things 
ought to be done. The new analytical power may well have more 
impact on future sampling and testing, and our basic under- 
standing of geotechnical processes, than is anticipated 


today. 
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APPENDIX A 


COMPARISION OF RESULTS OF FENA.3D WITH LEFEBVRE ET AL (1973) 


From the figures herein, it is seen that in general 
the results of FENA.3D agree well with those of the 
reference. Stresses plotted from FENA.3D are based on 
stresses in the elements adjacent to the longitudinal and 
transverse sections. The writer is unsure as to the plotting 
procedure used by Lefebvre et al. 

Particular attention is drawn to the stress and 
displacement results in the elements adjacent to the valley 
wall. Here the finite element formulation should lead to 
an artificially stiff structure. Lefebvre et al. discuss 
this situation and indicate that they discarded results 
from triangular elements at the valley wall for their 2D 
analyses. The improved accuracy of their 3D element may 
justify their accepting the results of tetrahedral and 
prismoidal elements at the walls, but experience with 
FENA.3D convinced the writer that boundary stresses must be 
accepted with caution. The shear strengths mobilized in 
these elements tended to be high and in the analyses of Mica 
some shear failure occurred. Unfortunately the doubts with 
which these results can be applied to the real dam prevent us 
from examining the development of stress transfer at the 


abutments in greater detail. 
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Discrepancies in the results can be attributed to 
aecombination of: 

(i) slightly different elements 

(ii) logic of the layered analysis 

(111) artificially stiff elements in contact with 
the valley wall. 

(iv) different program algorithms and computer 
systems. 

Results of the analyses indicate arching in the lower 
central region of the dam of about 28% for FENA.3D and 25% 
for the other, for the 1:1 valley case. Mica Dam has roughly 
2:1 valley walls for which the results of Lefebvre et al 
indicate about 18% arching. The homogeneous analyses of 
Mica Dam using FENA.3D yielded about 30% arching due to the 


3D valley geometry. 
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0.5 FENA.3D 
LEFEBVRE, ET.AL. 


VERTICAL SETTLEMENT 


0.4 FENA.3D 
LEFEBVRE ET.AL. 


HORIZONTAL MOVEMENT 


TRANSVERSE SECTION: HORIZONTAL AND VERTICAL MOVEMENTS 
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APPENDIX B 


DETAILS OF THE MICA DAM MESH 


Data for the program FENA.3D is prepared as per 
the users manual (Krishnayya, 1973b). The mesh described 
in Chapter IV is given in two forms: 

(1) the outlines of element boundaries for the 
different sections are given in Figs. B.1 to B.9. Solid 
lines are in the vertical plane of the section and dotted 
lines are not. Nodes are numbered consecutively through 
the set of planes so as to minimize the bandwidth of equations. 
Elements are numbered consectively within each layer. These 
sections can be used in conjunction with, 

(2) a full set of nodal data and element data as 
required for the program FENA.3D. The space on the figures 
is limited and to avoid confusion only essential numbers are 


given. This set of data follows the figures B.1 to B.9. 


ee 


(= 
wen ns 18 : 
ne a. mses oe 


tsdfizbeb dgom ent . (devel , 
sami02 ow? i aeviv ea vr 
3 36 eeiinbaued anomie to sontlsuo arid 4 
hrle® -@.4 63 1.@ .aeplht ai novip o3e anolscee . 
besteh bot aotsoee afd fo eanaig [soittev ofts a 
tipvouts ylovicuovede Serodmn 916 zabow 
ano) Jeups 30 Acbiwhnad oft sgiminim of as oe aensig 3a ean a : 
Sean? .isys! dees nidviw ylovitoeseqco baisdmpn 916 = 
<3 iw nolsoauinos at beet ed mso sz 
as 6fah dnometp bas e¢6& [ebon to tee Ilu3 s (S$) 
eoudpti odd ne sosge edT Of .Axaa msxpowW end 202 Bex! 
exe a'actewn igiineees ytdo notautces bievs e3 Bas bastahl 
.2.8 of 1.8 zeunpit erg swollo3 s1sb to 


ATEN) 


O71Z~ 


o0tZ 


0052 


99” 


INAWINGVY LHDIY “ONTdOTS :1 NOILOSS 1° FUNDIA 


o92 


902- 


9S 6) 
a / if Wa 
/ / ve 
mg of. A 
7 y, / 7 
oo / / pie eee 
ie en 
/ ee 


oob- 


009 - 


908 - 


at 

oe - 
a \ a a 7 
7 ae iy aa 


oa By ie 5 7 


mas 
ey hamec ee 
ae sua eae 
- 


| Nemes 
=a 
Sh ae 
ate 
1 ae = oe 


Sw He, 


O¢+ST * VIS 


*¢ NOILOES 7@°a HUNDIA 


149 


06+9T*° VIS *€ NOILOUS E€°RA ANNDIA 


EIChbE 8°? SRCLIOd 3+ ery’ TessG 


* 


Bieri ta SES 


- 


< 


a G 
ae eee te 


an 


-~. ; 


= a ocr Ss = 


141 


OO0+1Z°VLIS ‘7 NOILOES 7°a HANOIA 


Pro rer); ae1 oF 


fo o}=)) 099) 90%! Qo7) 099) 00g 


aes ieee: Sars mci Sat (eon ee Se Ga ee ee i 


210) oo coz ° OCe~ Cee 225 00g- ecg) — Q07}~ 


~ 


9orl- 


T 


1 


4 


- on marae 


kK ig = Gr ae 
7 


142 


OS+2¢°VIS *S NOILOgS 


G°@ Wana 


- 1 | | peas dt 1 ie a es 


ll 


~~ = a 


— _ i adie 


Oeecs AT :& WOTTIe 3c.) 2ab0re 


143 


OS4+G67°VLS ‘9 NOILORS 9°A WANDIA 


99 5 Oo>! Qo7Z! go9) 902 299 oob lord ovz~ Oey Q09- ocg- cool- ooz]j— O0%)— 099)~ 
4 


| 00S) | 


091% 


144 


0S766 Vise 2 NOTIOIS * 2° Aan 

ro 2° fo i a a i” a "oF ine | 

02g asg Qo Qoz : reley am Cop 999- o93- 999) — goTl- 20n1— 2991- 208/— = 
TOOL! = 
— 926) | 
be 4 
r— aoe a 
t— 99% = 
é Z 
= o9sZ 3 

i ! | Ie | a | za] = 1 


O¢4et. ate <\ MONTOSe 1,8 2MIS 


145 


O8+€E°VIS *8 NOILOUS 8° AUNDIA 


Q07 


-2- 


0903 — 


4 


akinesia ot acai Cea on 


Q090)— 


O9Z)- 


20) - 


09 /— 


oog)— 


a ig 


O°00Z~ 


ae INAWINGVY LAT “ONTdO1TS :*6 NOILOUS 6° FUNDA 
| 
| a es em Gar RY Wie Pa ae er ate ie 
a eo) oor aor o 202— Cora 209- 293- 099)— 90° Or! — oo9g)- 20g)- O00T~ 
; Ook) 
-— 906/ ee ale aise ene y Zz 
= oe 
oe PY; 
I— oolz ee 
L ye 
ame OS estes 
raaperA \Lz azr b5 
L 
| l fe gee ck, Seton old Ge Me ag al ee Sp a pn ee Bs 


MICAS DAMP STRESS 


206 254 
NO. KODE xX 
1 12 6=59 
2 12 35 
3 72> AGG 
4 q2 -320 
a 12 -260 
6 127 2-183 
7 120 “#8 
8 42 435 
9 42. 938 

10 42. 1246 
14 () 1565 
12 0 =6 
a3 0 54 

14 0 -521 
a6 b =280 
16 0 +128 
17 OG «= 86 
18 0 .=46 
39 a. Jd57 

20 0 284 

21 A 7 

22 12..=449 

23 12> —=205 

24 AZ =142 

25 12) «395 

26 12 9AGd 

27 To) G4 

28 1Z 196 

292 12 38A2 

30 12-76 

31 tant 3 2G 

32 1264209 

33 1é, 3445 
34 12 «270 

35 12 140 
36 S 2th 

37 io 2e4 

38 9 34 
39 G -338 

4o O  =235 

4 0 -140 

42 0 -104 

43 0 66 

4u 0 136 

45 % 265 

46 12> 1020 

47 12 -~749 

48 mW =3 75 

49 0 -206 

50 O = 145 

51 QO +483 

Sz ha 9930 

a Va VSi2 

54 12, 90 


7 
4 
1240 
1240 
1246 
1400 
1470 
1470 
14806 
1490 
1480 
1455 


1690 


254 


Z 
2500 
2500 
2500 
2305 
2320 
2320 
2320 
2320 
2876 
2400 
2500 
2500 
2500 
2295 
23.20 
2320 
2320 
23.26 
eae0 
2345 
2180 
2255 
2276 
e205 
2292 
2300 
2300 
2300 
2423 
2160 
2145 
Pea er 
Z2U55 
2455 


ANSLYSASa- 5 LIFTS 


U 


0 


147 


Sets gaet Oe> St Tom 
#2rs Ofat 00S- Sf SE -) 
a 


sek ed ek ek eed de ed wd 
ONMNNNNNNNNHNNNNNNN NW NH 


saath: sno 
OO OOOO WNnNnN OOS SG GOOD OoOODO OOO Os 


1690 
1690.2 


1690 
1690 
1960 
1675 
13830 
1840 
1845 
1800 
1770 
2060 
1985 
2620 
1925 
1890 
1920 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
21090 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
21006 
2100 
2100 
2100 
2100 
2100 
2100 
2190 
2100 
2100 
21090 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
2100 
22596 


148 


: inenenececenronccssobosness i 


rer erste iar ree ad | 
388i ya gpeusasipHiTSHET TEESE 7 | 
PESEEEEE rit: pris aucun sEaEESaS 


org ea 5 
0 6 
0 -568 
0 -3206 
0 -166 
0 -128 
0 -90 
O F490 
0 237 
0-1010 
0 -784 
0-414 
0 -232 
0 -170 
0-168 
O=Thay 
0 464 
0 652 

12-1330 

12-1113 
0 -625 
Q -3909 
0-211 
0-124 
0 284 
0 767 
0 944 

12 -770 

12 -468 

12 -298 

12 -150 
0 80 
0 438 
0 820 
0 = 375 

12 360 

12 900 

12 1330 
0 1335 

12 1700 
3 1700 
0 -106 
0 -50 
0 4 
0 -562 
0 -325 
0-170 
0 -134 
0 -96 
0 104 
0-234 
0-1010 
0 -786 
0 -445 
0 -236 
0-174 
0 -113 
0-120 


2250 
2ee0 
2250 
2250 
Deo 
2250 
2250 
2250 
2250 
2250 
2250 
2250 
2e2au 
2250 
Ze2aQ 
2250 
2250 
2250 
2250 
ees) 
2230 
2250 
Page 9) 
Ae} 0) 
2250 
Pig ace ih 
2250 
2250 
2250 
2250 
aay 
2250 
2250 
2250 
2250 
2200 
2250 
wand 
2250 
Z250 
a dre 5) 
2550 
RooU 
2550 
2550 
2550 
2550 
2550 
2550 
2550 
25500 
2550 
25590 
A590 
2530 
2200 
2550 
2550 
ie be 
2550 


2500 
2500 
Ze 
2320 
2320 
2520 
Zoe 
PLENER® 
2345 
2095 
2114 
2140 
2155 
2195 
2155 
2150 
2143 
2738 
1350) 
1964 
1996 
1990 
1990 
1990 
T7490 
1990 
19:90 
1885 
1818 
1746 
sas) 
1830 
1810 
1810 
1885 
1700 
peta) 
1720 
1885 
1740 
1860 
2500 
2500 
2500 
2296 
2320 
2320 
2320 
2320 
2320 
2345 
2096 
2114 
2139 
2155 
21905 
2155 
Z1sg 
2Tti2 
are 


149 


7) 
+ 


1 fe: ew 


"a. - 
> 7 


SRSA RARER ES A vag ESR Ws 


SE5S% 

4436 < 

—-4 = 7 
r—% 2 
oe ee 
el ah ae oe 


oyyunasae 


hae tk ea 


seetaateapitstGHe 


aouweo 
23 i3eas 
533% 


=. — 

- eal 
Biebited 
a 


Beans 
: ia 


a 
= * 


> ¢ \~ 
= be < 


0-12506 
0-1063 
CG -674 
Ge -303 
o -24'5 
0» =136 
f= §23G 
12 “7 8G 
12 %6!2 
12-1690 
A= 1538 
12-1928 
2) =3 90 
0 -686 
12-7000 
12 Too 
eo Sac 
42 —26) 
Va He 
42-305 
12. 340 
12 868 
12 1020 
12 1320 
12° 7366 
12 “259 
$2, 1826 
10 1650 
12 1650 
Gr h=95 
60 -40 
0 15 
o-=5 39 
o~=30S 
oe =150 
oo =Age 
G =t5 
0 116 
0 243 
60-1000 
a =775 
gy =8 55 
225 
0 -164 
v= 3:02 
12> 73a 
te) ARF 
te. 539 
G- 1240 
0-1048 
2-670 
12 -400 
12 -240 
Jee =a 1e 
ta 290 
12 490 
12 646 
o- 1650 
0- 1480 
O0-1114 


D0 


0 -989 
9 -692 
12 -530 
12-1659 
12-1480 
12-1210 
12-1000 
12 -695 
o =70 
0 -10 
0 50 
12 -646 
12-370 
12 -132 
12 -86 
12 -42 
ta Sg 
12 300 
0 -880 
12 =—353 
12 -448 
12 -250 
12 ~162 
12 -88 
90-1214 
12 -992 
12 -660 
3-1700 
00-1510 
12-1186 
12 -990 
12 -690 
12-1700 
12-1506 
12 -130 
12 22 
12 160 
12 -690 
12 -955 
12-1199 
412-1710 
12-1470 
KODE Nl 


1@3 


WWNHNNDN WN WW WNNDN DN ND LO 
anh 
i) 
Ne) 


bbe l 


MATNO. 


—_ 


ee ee eee eee ee eee eee a ee ee ee ee ee ae ee 


; 


faa > 


4 _ 
a. 
rs 


nbehysongese BGs 


ta. | 


reapTE 


e228 
ca) 


bs 


NNNNNNNNNMNWWWWNYNYWWNWWWDNYWNHYWWHWWWNWNWWNYNYDNNNYNNNHNNNMNNWNHY WYN MYM WNH Ww 


SSP SS PSS SEMA Tenuta Mle SS PS) ICS ESS Je ES IG aS aot ne ete Ore ee fe Salerro arene | om em ee CO eres ROU al fo neta ye eal) men ean re BeNOR el ee ly SNe Nay ay Gon son 


132 


7 ES a 


a et es Ee | Uwvoua a 4— e745 w=" = 3 <3 2 


ee ee, oe | a) 


Maa Aas Oa AI Aaa 


guaesaaaeecetsccseszaisscaasteoeesecsesRaness 


- 
i — 


MPN NM FWWHONNNHN EF EWWNHNNDYN EF WWNNNDN SWWNNMNHDHNGS EF ES EWWWWNHNNDN ES SE WWH WD I~ 


153 


ie 


AEM MAME MM EST EMH eR PAH Hh 6 ERMA 


oe ae ie ad alt cl at oa ro eee 


4 
ae 
ae 
ee 
s? 
Be 
Te 
ap 
ett 
wer 
af 
oe 
ver. 
BEr 
oe 
ovr 
Ter 
atr 
err 
oar 
rar 
cer 
one 


n 
3 
aa 


SERESSIEISIS SSSR ASASIALATS 


ae 


~~ 


is 


bade 


_? - 


{= 
_ 


AOR Re mmm mm or 


} 


38 


144 
145 
146 


194 
195 
145 
a7 


a geal SET Seung SEP E Ne CEE meee crves CEerw reer erer ere tree 


wen eak ln ett eh aa eS I sh. Bo yh Nd tio ty a hg nd I a, ta ta) AL Rab Ea gee Dada) Ba ata N UI Nee RO ORRSe aE ERS 


AnaWAte eS ee MAME 


ote Te = ris 6? 


ces vee 


mins eSrhnmnwnear pr A 4 As 


) sa tf 

we te 
ts 4 

st 20 arti 

ow #8 e 
ts ar 

ke er | 6 

e¢ re 

ee oy ph 

£ 

te c 68 
ri TY v te 
ed rh 

es @F  JEY 6&8 
ty or 486eF te 
en = et ep 
Pa t “a pe 
at atr F977 et 


a= 


et 
et ve 
oe eT 
fa 08 
arr Tre 
err arr 


2 
(26 “2€ a 
ne 62 Ct a 
te c8 38 © 
ge $4 ‘€¢ Sa € 
ea oa ‘ea Ss 
oc B68 6 eB o et 
Tt? Te a6 re 
se 68 a «f2 93> 
ea Be. *S? cf 
oe “ee  e@8 e@2 "ft 
EA ast er’) s6 oe 
os esr “2s! ©&e is 
as rar | oer 966 . 
an ‘esr War eme . €& 
te ecr r - et 
an 40k est eee  § 
on 4 6ver oer 88 4 E 
© er ote 
tr er => 
act Bar oer 
gar asr’ & 
Rar die 6 
iFSe 
gr. 
Lf & 
rer € 


sean 


es 
ae 
" a + 


m annmnn Talos 


NNNNNNNNWWNHYNNWWNYNNNWWNHNYNNWWNHYNNNWWHNYNYNWWNNNNWWNMYNNYNYNY WWNND 


247 


NONMNWWNHONYNNWWNYONYNNWWNYNYNNWWDNNYN WW DONDNWWNHDDNWWNNDNNWWPNDOND &WWNDN 


Beye 


VEU UWE EEUU We UOTE WeMouuuwunbaNaamny ‘ 
reeneenstnenres tine 
geaseStaseescsss pained seinneee T17 
weesesraitssises: 


; % it Se pape te fo Ps 
S53; Seisowserassnnsses 


ae 


BeSgsSe  ouutzas ceweceneskers 


Aa A GAN GA ANN ak 8 A th el ss St ta ae at se tt te AA et a ee Ho A es COE A 


APPENDIX C 
DETAILS OF MV GAUGES 


On succeeding pages are some instrument details 
supplied by CASECO. The first shows the installation and 
layout of the telescoping-sectioned MV gauges themselves 
(Fig. C.1). The second shows the layout of the latch-cone 
device used in conjunction with a survey tape to measure 
casing profiles for settlement and incremental strain 
data. A similar device was used with a hauling cable in 
horizontal (MH) movement gauges, which were themselves 
Similar in most respects to the MV gauges. 

Figure C.2 is a photograph of the latch-cone device. 
In the closed ‘position, at right it latches under each 
section of casing. When it reaches the bottom of the gauge, 
it is dropped to jerk out the springs, as shown on left, 


that enable it to be drawn up out of the gauge casing. 
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APPENDIX D 
PLOTS OF STRESS DISTRIBUTIONS FROM FENA.3D 


Herein laretavseries sofuplotsrosistresses cas 


described in Section 6.3. 


Transverse Section, Station 22+50 

Plots show major, intermediate, and minor principal 
stresses and maximum shear stress for the following cases: 

(a) 3D linear analysis, 7 materials, 1 lift, bedrock 
ae ed 1584 

(b) 2D linear analysis, plane strain, 7 materials, 
5 lifts, bedrock fixed. 

(c) 3D linear analysis, 7 materials, 5 lifts, bedrock 
settlements included. 

(d) 3D multilinear analysis, 7 materials, 5 lifts, 
bedrock settlements included. 
As mentioned previously, the stresses shown are those computed 


at the element centroids between Stations 22+50 and 25+50. 


Longitudinal Core Section 


Plots show major, intermediate, and minor principal 
stresses and maximum shear stress for cases (a), (c), and (d) 
detailed above. As mentioned previously, the stresses shown 
are the mid-core average of values computed for the upstream 


and downstream core elements ata certain section. 
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The reader is asked to interpret these results in 
a suitable fashion. Clearly values for the incremental 
analyses do not differ markedly but the 1 step "gravity 
switch-on" stresses are not so agreeable. In general, 
major principal stress tends to align with vertical stress, 
intermediate principal stress with cross-valley (Y direction) 
stress and minor principal stress with transverse (X 
direction) stress respectively, although the trajectories 
curve in an appropriate manner near boundaries (e.g. 
Covarrubias, 1969). 

The writer concludes that the stresses from the 
linear analysis including pedvock set lenente are the most 
representative since the best agreement with field deforma- 
tions was achieved in this computation. The multilinear 
analysis has too soft a response in the 1969 core material 


(below el. 1990 f£t.). 
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TRANSVERSE SECTION: MAJOR PRINCIPAL STRESS; 3D, MULTILINEAR 


a 


AABATIZEIOM .CC ;2SNTS QAMPIMISY MOlLan :MOT Ose FeRTVSCART 


16Z 


PO 


NOOUCHE+UVANIT “dE ‘SSHYLS TydIONIYa YOrVW *NOILOGS ASYAASNVUL 


MI SANO LUO) 


L638 


‘5 ul SANJA %) 


NIVULS ANVId ‘dz 


*SSHULS TVdIONIYd YOrVW ‘NOTLOAS ASUAASNVUL 


WIARTS MMAI 10. ~eehsre IAT AST kotat 


164 


IdIT 1 ‘UVENIT ‘d€ ‘SSHYLS TVdIONIYd WOFVW *NOILOAS ASUTASNVUL 


{sy WL Sunojvo) 


_ 


Tihs I AAW OE [ERNST UaetoOWIAT nota -WORTDaS Sapavemeat 


165 


UVANITILIOW “d€ ‘SSHMLS TvdIONIUd 


‘$s UI $400) U07) 


ALVIGHWHHINI *NOILOSS ASYHASNVAL 


- 


—. 
= 


AABRELITIUM (GE ;22UeT2 SATIS? STAIGINATTAI -WOLTOUE S2mIVERART 


166 


MOOUGAAHAVANIT ‘GE ‘SSHULS TVdIONIUd ALVIGHWAALNI :‘NOILOAS ASUHASNVUL 


}s> uw S500}U0) 


\ 

\ 

' 

| 

‘ 
2 


a 


~ 
es ~ 


Low 


NIVaLS ANVId ‘dz ‘SssH¥uLs TvdIONTUa ALVIGHWHAINI -NOILOGS ASUAASNVUL 


ys VI SuNOpVO) 


WEARTS SUAIS 


® 


1 ,ateeae tA) DWI A 


iTalGenmiTn 


‘“OLYaee Fens eRAet 


168 


TAIT I ‘UVENIT ‘GE *‘SSHULS TVdIONIUd ALVIGAWHAINI *NOILOAS ASYMASNVAL 


‘Js UV! S4nojzvoy 


_sipeme nee EPoR Ow: IWLEHROIVAR BEINGIEYY BieReet 3° rraeye’ J Tres 


169 


UVANITILIAW “d€ ‘SSHYLS TVdIONIUd YONIW :NOILOAS ASMAASNVUL 


‘$5 ul Ssnojuo*) 


>, 


_ _ 
: 
= | 
Ee - 
_ 


AML TAM HE -agTNTe cATTOMTAT ROWE MOL eR 
_ 
SenavaMaet 


170 


MOOUCHAHAVANIT ‘GE SSSHYLS TVdIONIUd YONIW 


+ vw oor) 


*NOILOAS ASYHASNVUL 


Begin 


NIVULS ANVId ‘dZ ‘SSHULS TVdIONIYd YONIW :NOILOUS ASYHASNVAL 


ys SINOJMO’) 


WIA aWAIY US )32997> sAqnD4ras HanTM -nolTode SOaaveWAAT 


By) 


LdIT T “YVaNTT ‘ae 


‘Jsy WL suoojue) 


¢ 


*SSHAULS TVdIONIYd YONIW *NOILOHS ASUMASNVUL 


yarpes’ 


SAARP2E SECLTOW- wieoe bite, vT 


7 


Las 


UVANITIIIAW “GE ‘SSHYLS UVAHS WOAWIXWN *NOILOGS ASUAASNVUL 


J$> Ww suNg|Uo") 


AAT TUOM GE ;eeANTS ANQUe MMEXAM “Mart ake BaesveMert 


174 


MOOUGHATAVANIT ‘dE 


fs Uy Sue o|U07) 


‘SSHULS UVAHS WAWIXVW *NOILONS ASYAASNVUL 


he 


hae 


uw 


NIVULS ANVId ‘GZ ‘SSHULS UVAHS WAWIXVW ‘NOILOWS ASUAASNVUL 


s3nouo’) 


176 


LIT 1 


3s) uv) SANOL UO") 


¢ 


UVANIT ‘d€ SSSHULS UVAHS WOAWIXVW :NOILOES ASUSASNVUL 


1g 


UVANITILION “d€ !SSHUYLS TVdIONIYd YOFVN *NOILOAS TYOOD *LIONOT 


‘\s> ul SANO JU’) 


x) 
= 
ro) 


OI1+8€ 


178 


O1+dl 


MOOUGAGHAVANIT “dE ‘SSHUMLS TVdIONIYd YOFVW *NOILOAS AYOD °LIONOT 


Js Yl sunoqjuo’) 


01 +95 


17.9 


Ol+el 


IIIT T “UVANIT “dE ‘SSHYLS TVdIONIYd YOFWW ‘NOILOTS FAOO *LIONOT 


JSy YU! Sunoyu07) 


01+8¢ 


fuetE” COKE ECR NOW! weigr Wuernctayr eteeee!. so" reso” 7 ORS 


180 


O14! 


MVANITILION ‘de 


*NOILOES FYOD *LIDNOT 


*SSHULS IVdIONIYd ALVIGHWAALNI 


‘sy UI ssn} ’) 


— 
— — —— 
—— oe eee 


—_ 1~—. — 
eee 


01+8E 


— 


CSR SEOLYOK! LAECEWMSIOLE be betSeyT aLinee* a’ WHITE laeve 


131 


Oj +2) 


MOOUGHE+AVANIT ‘GE *SSHULS TVdIONIUd ALVIGAWAAINI 


>NOILOES HYOO °*LIONOT 


‘\sy ul s4nojyo~) 


O1+8E 


182 


Ol+el 


LALIT T 


— ee 


“YVAENIT ‘ae 


*SSHULS TVdIONIYd ALVIGAWYAINI *NOILOAS FAOO *“LIONOT 


JS) UL sungyuos) 


01l+8E 


pens) 


MVANITILION ‘G€ ‘SSHUYLS IVdIONIYd WONIW *NOILOAS AYOOD *LIONOT 


Jsy Vl sunojuo ) 


Ol+2! 


OI+8€E 


WACULETLUM (GE eNORVe SES 14IAT AORIM -HON TERS F800 JTC 


184 


MOOUGHAHUVANIT ‘GE ‘SSHULS TVdIONIUd WONIW :NOILOS FHOO * LIONOT 


ot oe ul sunojuo’) 


ro 


a 
— 
—_—— 
ee eee 


Ol+al 
o1+8E 


TOuetL” Corp akcRiou: arioscvivciiwr cayEsa’ 20‘ PieReeRuENEK: 


185 


LATT 1 


‘MVANIT ‘dE 


*SSHULS TVdIONIYd YONIW :NOILOAS AYOD *LIONOT 


‘|sy u! SUNOIUO *) 


01 +8¢E 


* Cope Seeryout wrreor wemibawr ereetat do" revere’ 7 Prag 


Fruye tt 


186 


UVANITILION “dE ‘SSAYLS UVAHS WAWIXVN ‘*NOILOWS FYOD °“LIONOT 


1s) ul SANO JUO’) 


O1+ZI 
01 +8&F 


187 


O1+d2l 


MOOUCAA+UVANIT ‘GE ‘SSHULS UVAHS WOWIXVW ‘NOILOAS AYOO °LIONOT 


JSD}. Ul SuMejuO>) 


o1+8€ 


FOUMNEAR BART] . oe 2 52RATE DAS MOMIZAM : WOSTGRS WROD , Tea 


188 


Ol+2\ 


LATT T 


“AVANIT ‘GE *SSHMLS UVAHS WOAWIXVW :*NOTLOHS AYOO ° LIONOT 


ae 


ul S4N0jUO™) 


ol+8E 


APPENDIX E 
COSTS OF THE FENA.3D ANALYSES 


Apart from the previously mentioned data preparation 
costs, the chief costs in computation are the generation 
of structure stiffness and the solution of the large set 
of simultaneous linear equations to obtain nodal displace- 
ments. Structure stiffness is composed piecewise of element 
stiffness, thus the cost of this factor is proportional to 
element type, integration procedure, and the total number 
of elements treated, The solution time of the equations 
depends upon the semi-bandwidth of the set, which is a 
property of the mesh of elements. Other costs are 
marginal, and for evaluation purposes the cost of an analysis 
is a function of total elements generated and equation semi- 
bandwidth. 

Krishnayya (1973b) presented a correlation of semi- 
bandwidth and CPU time for solution of the equation set. 
This is shown in Fig. E.1. It must be remembered that, as 
used, the program FENA.3D depends heavily for its economy 
upon MTS system subroutines set up for the IBM 360/67 at the 
University of Alberta. The writer found a more useful 
prediction was the total CPU time for a run, correlated with 
the product of semi-bandwidth and total elements. The 


relative cost of CPU storage increases with the overall core 
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area required for a problem too. Some costs of various runs 
for the Mica Dam analyses are given in Table E.1. 

Finally, an information sheet conforming to the 
requirements for the W. E. S. Symposium (Desai, 1972) is 


given. 
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FIGURE E.,1 APPROXIMATE CPU TIME FOR SOLUTION OF EQUATIONS 
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COMPUTER CODE (PROGRAM) QUESTIONAIRE HSE) 


Name of Code, if any: FENA . 3D 
Purpose of the Code: To compute stresses and strains using linear 
or non-linear stress-strain relationships. 


Developed by (including year of development): Krishnayya, A.V.G.(1973) 


Modified by or revised by (including year of modification or revision): 
Simmons (1974) for multilinear analysis and tapewriting options. 


Code available from: Dept. Givil ine., Univ. of Alberta, Edmonton. 


. Brief description of the problems solved by using the Code: 


Three dimensional constructional and cracking analysis of earth 
dams; Building settlements, Excavation movements. 
Other types of problems for which suitable: Any three dimensional 


problem with linear or non-linear stress-strain relationships. 


Type of finite element used: Isoparametric hexahedron with 8 nodes. 
Type of nonlinear technique used: Piecewise incremental procedure, 
Type of displacement function (or field variable function) used (e.g. 
linear, quadratic, etc.): Linedr. 

Type of equation solver used: Gaussian reduction and back-substitution 


lig, lop ilfeyellkss\ 


Language used: FORTRAN IV. 

Machine on which presently adopted: IBM 360/67 with MIS operating 
system. 

Capacity in terms of number of nodes and/or elements: Maximum elements 
350, Maximum nodes 450, Maximum semi-bandwidth 250. 


Is documentation such as Guide for Data Input, etc. available? Yes. 


10. Other comments such as number of cards in the code, computer time 


required for a typical problem, limitations if any: 1545 cards. See 
Table E.1 for guide to times. Uses 8 MIS system subroutines to 
achieve economy. Requires two temporary disk files, optionally 
also one tape-drive, Data preparation time and checking procedures 


may be measured in months. 
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